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ABSTRACT 

In order to determine the physical properties of the hottest and most luminous 
stars, and understand how these properties change as a function of metallicity, we 
have analyzed HST/UV and high S/N optical spectra of an additional 20 Magel¬ 
lanic Cloud stars, doubling the sample presented in the hrst paper in this series. 

Our analysis uses NLTE line-blanketed models that include spherical extension 
and the hydrodynamics of the stellar wind. In addition, our dataset includes 
FUSE observations of OVl and HST near-UV He 1 and He 11 lines to test for 
consistency of our derived stellar properties for a few stars. The results from the 
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^Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer, operated 
for NASA by John Hopkins University under NASA contract NAS5-32985. These observations are associated 
with program C002. 

^ Visiting astronomer, Cerro Tololo Inter-American Observatory (CTIO), a division of the National 
Optical Astronomy Observatory, which is operated by the AURA, Inc., under cooperative agreement with 
the National Science Foundation. 



complete sample are as follows: (1) We present an effective temperature scale 
for O stars as a function of metallicity. We find that the SMC 03-7 dwarfs are 
4000 K hotter than Galactic stars of the same spectral type. The difference is 
in the sense expected due to the decreased signihcance of line-blanketing and 
wind-blanketing at the lower metallicities that characterize the SMC. The tem¬ 
perature difference between the SMC and Milky Way O dwarfs decreases with 
decreasing temperature, becoming negligible by spectral type BO, in accord with 
the decreased effects of stellar winds at lower temperatures and luminosities. The 
temperatures of the LMC stars appear to be intermediate between that of the 
Milky Way and SMC, as expected based on their metallicities. Supergiants show 
a similar effect, but are roughly 3000-4000 K cooler than dwarfs for early O stars, 
also with a negligible difference by BO. The giants appear to have the same ef¬ 
fective temperature scale as dwarfs, consistent with there being little difference 
in the surface gravities. When we compare our scale to other recent modeling 
efforts, we hnd good agreement with some CMFGEN results, while other CMF- 
GEN studies are discordant, although there are few individual stars in common. 
WM-Basic modeling by others have resulted in signihcantly cooler effective tem¬ 
peratures than what we hnd, as does the recent TLUSTY/CMFGEN study of 
stars in the NGC 346 cluster, but our results lead to a far more coeval placement 
of stars in the H-R diagram for this cluster. (2) We hnd that the wind momen¬ 
tum of these stars scale with luminosity and metallicity in the ways predicted 
by radiatively-driven wind theory, supporting the use of photospheric analyses of 
hot luminous stars as a distance indicator for galaxies with resolved massive star 
populations. (3) A comparison of the spectroscopic masses with those derived 
from stellar evolutionary theory shows relatively good agreement for stars with 
ehective temperatures below 45000 K; however, stars with higher temperatures 
all show a signihcant mass discrepancy, with the spectroscopic masses a factor 
of 2 or more smaller than the evolutionary masses. This problem may in part 
be due to unrecognized binaries in our sample, but the result suggests a possi¬ 
ble systematic problem with the surface gravities or stellar radii derived from 
our models. (4) Our sample contains a large number of stars of the earliest O- 
types, including those of the newly proposed 02 subtype. We provide the hrst 
quantitative descriptions of their dehning spectral characteristics and investigate 
whether the new types are a legitimate extension of the effective temperature 
sequence. We hnd that the NIII/NIV emission line ratio used to dehne the new 
classes does not, by itself, serve as an effective temperature indicator within a 
given luminosity class: there are 03.5 V stars which are as hot or hotter than 
02 V stars. However, the He I/He H ratio does not fare much better for stars 
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this hot, as we hnd that He I A4471/ He H A4542, usually taken primarily as a 
temperature indicator, becomes sensitive to both the mass-loss rate and surface 
gravities for the hottest stars. This emphasizes the need to rely upon all of the 
spectroscopic diagnostic lines, and not simply N HI/N IV or even He I/He H, for 
these extreme objects. (5) The two stars with the most discordant radial veloci¬ 
ties in our sample happen to be 03 “held stars”, i.e., found from the nearest OB 
associations. This provides the hrst compelling observational evidence as to the 
origin of the held O stars in the Magellanic Clouds; i.e., that these are classic 
runaway OB stars, ejected from their birth places. 

Subject headings: stars: early-type, stars: atmospheres, stars: fundamental pa¬ 
rameters, stars: mass loss 


1. Introduction 

An accurate knowledge of the stellar ehective temperature scale provides the means 
of converting the observed properties of a star into its physical properties. The issue is 
particularly important—and challenging—for the massive 0-type stars, as the bolometric 
corrections (BCs) are a very steep function of the assumed ehective temperature (Tefj) at 
the extremely high Te^’s that characterize these objects. For these hot stars, a 10% error 
in Teff results in in an error of 30% in the derived bolometric luminosity of a star (starting 
from the absolute visual magnitude; see Massey 1998), compromising any attempt to use 
stellar evolutionary tracks to determine distances, initial mass functions (IMFs), and ages of 
clusters (see, for example, Massey 1998, Slesnick et ah 2002). In addition, a 10% uncertainty 
in Teff results in a factor of 2 or more uncertainty in the Lyman continuum hux, ahecting our 
understanding of the ionization balance of H H regions and the porosity of the interstellar 
medium in general (see, for example, Oey & Kennicutt 1997, Oey 2004). 

At the same time, the determination of the ehective temperature scale of 0-type stars 
is complicated by the fact that the stellar atmospheres of these stars are physically complex 
(see, for example, Kudritzki 1998). The strong spectral features are formed under non¬ 
local thermodynamic equilibrium (NLTE) conditions (Auer & Mihalas 1972), while stellar 
winds provide a signihcant source of heating of the photosphere through the backscattering 
of radiation (Hummer 1982; Abbott & Hummer 1985). Modern stellar atmosphere models 
for 0-type stars include NLTE, spherical extension, the hydrodynamics of the the stellar 
winds, and the ehects of line blanketing; see, for example, Pauldrach et al. (2001), Hillier et 
ah (2003), and Puls et al. (2003, 2005) for descriptions of the stellar atmosphere programs 
WM-BASIC, CMFGEN, and FASTWIND, respectively. 
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The most recent improvement in these codes is an improved treatment of line blanketing. 
The results of htting these improved models to stars have suggested that the previous effective 
temperature scale for Galactic O stars may be too hot by as much as ~ 20% (see, for example, 
Martins et ah 2002, Bianchi & Garcia 2002, Garcia & Bianchi 2004, and Repolust et ah 2004), 
while the results for the Magellanic Glouds are mixed (compare Growther et ah 2002, Bouret 
et ah 2003, Hillier et ah 2003, and Martins et ah 2004). 

We are engaged in a new determination of the effective temperature scale of 0-type 
stars, with an emphasis on how the conversion from spectral type to Tefj depends upon the 
initial composition of the gas out of which the star was formed; i.e., how do the physical 
properties of an 07.5 V star differ in the SMG, LMG, Milky Way, and the Andromeda 
Galaxy, systems which span a range of 10 in metallicity. In a previous paper (Massey et 
ah 2004, hereafter Paper I) we analyzed the spectra of a sample of 20 SMG and LMG O- 
type stars using FASTWIND. The same code has recently been used to analyze a sample 
of Galactic 0-type stars (Repolust et al. 2004), obtained at a similar signal-to-noise ratio 
(S/N) and spectral resolution (Herrero et al. 1992). Our study suggested that the Magellanic 
Gloud sample is 3,000-4,000 K hotter than their Galactic counterparts for the early through 
mid-Os, although the sample size precluded a more dehnitive statement. 

Here we extend these studies to an additional sample of 20 Magellanic Gloud 0-type 
stars, and consider the results from the complete sample of 40 stars. In Sec. 2 we describe 
the space- and ground-based data of the new sample, along with our reduction procedures. 
In Sec. 3 we provide the model hts to these spectra, determining physical parameters, and 
describe the spectrum of each star in turn (Sec. 3.1). Our data includes very high S/N 
spectra of stars of the earliest types, including the newly proposed “02” type (Walborn 
et al. 2002), and in Sec. 4 we provide the first quantitative classification criteria for such 
stars. Our analysis also allows us to investigate whether the physical properties of such stars 
are well correlated with the new spectral subtypes (i.e., are 02 V stars necessarily hotter 
than stars classihed as 03 V or 03.5 V?). Our analysis of the complete sample of 40 stars 
then allows us to compare the physical properties of such stars as a function of metallicity 
(Sec. 5), including deriving effective temperature scales for the SMG, LMG, and Milky Way 
(Sec. 5.1) O stars, comparing the Wind-Momentum Luminosity relationship (Kudritzki et 
al. 1995) to that expected from radiatively-driven wind theory (Sec. 5.2), and investigating 
whether or not the derived spectroscopic masses agree well with the masses obtained from 
stellar evolutionary models (Sec. 5.3). We summarize our conclusions and lay out the need 
for future work in Sec. 6. 

We emphasize here that we realize that the model atmospheres are constantly improving: 
that just as the models today are considerably more sophisticated than the hrst NLTE models 
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of Auer & Mihalas (1970) on which the hrst modern effective temperature scale of O stars 
were derived (Conti & Alschuler 1971), so will the models 30 years hence will include better 
physics. Furthermore, even the models of today that purport to include similar physics do 
not necessarily yield the same results. So, we accept that ours will not be the hnal word on 
the absolute effective temperature scale. However, by using the same models and techniques 
on a large sample of objects in the SMC and LMC that have been used on Galactic stars, we 
can make the first cut at an effective temperature scale of hot stars that includes the effects 
of metallicity. We believe that the relative differences of the scales we derive are accurate 
and signihcant. 


2. Observations and Reductions 

The stars are listed in Table 1. As discussed in Paper I, good photometry provides the 
key to an accurate value of the absolute visual magnitude My, needed in order to constrain 
the stellar radius. Unless otherwise noted in the table, the photometry comes from the CCD 
UBV measurements by Massey (2002) or Massey et al. (2000), except for the R136 stars, 
which comes from the WFPC2 photometry of Hunter et al. (1997). Conversion to My for 
the non-R136 stars was done by using the color excesses determined from the spectral types 
(see Paper I), adopting distance moduli of 18.9 and 18.5 for the SMC and LMC respectively, 
following Westerlund (1997) and van den Bergh (2000). The My values for the R136 stars 
were derived in a similar manner by Massey & Hunter (1998). 

The spectroscopic observations are summarized in Tables 2 and 3. In general, our 
modeling efforts require data in the UV in order to determine the terminal velocity v^o of 
the stellar wind. These data we have obtained with Space Telescope Imaging Spectrograph 
(STIS) on HST using the FUV-MAMA detector. Data in the blue-optical region is used to 
constrain both the surface gravity g (from the wings of the Balmer absorption lines) and the 
effective temperatures Tgfj (from the He I and He H absorption lines). Since the detection 
of very weak (20-50mA) He I A4471 is needed for our modeling of the hottest stars in our 
sample, we chose to obtain long-slit data at high S/N. The long-slit observatioins allow the 
subtraction of nebular emission in this line, and we aimed for for a S/N that would allow 
accurate measurement and detection of this important He I line. For all but the R136 stars, 
we obtained 1.4A resolution data with the CTIO 4-m telescope and RC spectrograph, and 
a S/N of 400-500 was obtained with a great deal of care to flat-helding (see Paper I). The 
choice of spectral resolution was a reasonable match to the expected line widths. Spectra 
of the R136 stars come from CCD observations with RFF/STIS, obtained (in general) at a 
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S/N of 100 per 0.4A spectral resolution element, providing a sensitivity about half as great^. 
Finally, data at Ha is needed to determine the mass-loss rate M. For this we relied both 
upon ground-based CTIO and/or STIS observations. Full details can be found in Paper I. 
We note that in general our stars were isolated compared to the slit width (Table 2. Although 
one cannot in general exclude the possibility of an unresolved companion on the slit, these 
can be revealed by the sort of detailed analysis we undertake here. 

In addition, we also obtained data in two other wavelength regions for a few of our 
stars to provide some independent check on the modeling—how well do the data in these 
regions, not used for the model hts, agree with our results? One of these regions is the 
3020-3300A “near-UV” (NUV) region, obtained with HST/STIS-CCD in order to observe 
the He I A3187 and He H A3203 lines. We found relatively good agreement in Paper I 
for the two stars that were observed in the region. Of more critical interest is the far- 
ultraviolet (FUV) region. We obtained data in the wavelength range 905A-1187A with 
the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite (Moos et al. 2000). This spectral 
region contains the high excitation O VI A A1032,1038 doublet, which provides a useful check 
on the determination of v^o (see Figs 5 and 6 of Taresch et al. 1997). The large aperture 
(30” X 30”; see Table 2) necessitated the selection of isolated targets. 

The reduction procedure for both the HST and CTIO data are given in Paper I, and 
are not repeated here other than to note that the CCD spectra (both STIS and CTIO) were 
reduced with IRAF^ using the optimal extraction algorithms. For the i/FT/STIS CCD data, 
this led to a substantial improvement in S/N over that delivered by the pipeline (see Fig. 1 of 
Paper 1). For the STIS UV MAMA observations we used the standard pipeline reductions, 
as we found that reanalysis of these data made no improvement to that of the pipeline. 

The FUSE spacecraft and its instrumentation are described by Moos et al. (2000) and 
Sahnow et al. (2000). The FUSE spectra were extracted from the two-dimensional raw spec¬ 
tral images using CalFUSE version 2.1.6, which uses an optimal extraction algorithm de¬ 
veloped by S. Lacour (see http://fuse.pha.jhu.edu/analysis/calfuse.html). The observations 


^Others recently have made different choices in obtaining optical data for their modeling of MC O stars: 
Bouret et al. (2003) and Heap et al. (2005) rely upon the echelle spectra described by Walborn et al. (2000), 
obtained at higher dispersion than ours, but with a short slit that precluded easy sky (nebular) subtraction. 
The VLT-Flames survey by Evans et al. (2005) uses fibers on an 8.2-m telescope to sample many stars, with 
only modest compromise for nebular subtraction. We consider these approaches complementary to ours, 
each with its advantages and disadvantages. 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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consisted of from 2 to 6 individual sub-exposures, which were extracted, calibrated in flux, 
aligned in wavelength, and, if required, scaled in signal level to create a co-added spectrum of 
each star. Three of our stars (Sk—67°22, AV 476, and AV 177) were observed in July or Au¬ 
gust of 2002, following the discovery by the FUSE project of an error in the flux calibration 
procedure used by CalFUSE (see http://fuse.pha.jhu.edu/analysis/fcaLbug/fcaLbug.html). 
We have conhrmed that the effects of the reported error are entirely negligible (< 2%) for 
the three spectra at the wavelengths of the O VI lines, which are of primary interest to us 
here. 


These FUSE data are obtained in 8 channels (LiFla, LiFlb, LiF2a, LiF2b, SiCla, SiClb, 
SiC2a, and SiC2b). For four of the stars, the agreement between the channels was good in 
the regions of overlap, and we averaged the data to produce a higher S/N spectrum, binned 
at 0.1 A. For one of the stars, AV 177, the agreement between the 8 channels of FUSE was 
poor, with only two of the channels (LiFla and LiFlb) showing the expected flux levels and 
a normal spectrum. The observation was repeated at our request, but the new data suffered 
from an even worse problem involving data drop out, and we obtained useful data only from 
Lila and LiF2b. We combined these four observations for our treatment of AV 177. 


3. Analysis 

Following the procedures of Paper I, we first determined the star’s terminal velocity v^o 
from radiative transfer fits of the P Cygni prohle of the C IV A1550 doublet, and (where 
possible) Si IV A1400 and N V A1240, although the Si IV line is usually too weak, and the 
N V line is often contaminated by strong interstellar Lyo; absorption. The htting is done 
based upon the SEI method (see Earners et al. 1987, Haser 1995, and Haser et ah 1995). We 
list the resultant values in Table 4. We estimate our uncertainty in this procedure as 50-100 
km s“^, with larger uncertainties (up to 200 km s“^) indicated by Two stars had even 
less certain measurements, and we estimate the uncertainties as 300-500 km s“^. 

In running the FASTWIND models, we adopted these values of v^o and began with 
the assumption that (3 (a measure of the steepness of the velocity law) has a value of 0.8 
following Puls et al. (1996). In addition, we began by adopting a He/H number ratio of 0.1, 
and adjusted it if needed. The absolute visual magnitudes Mv' were drawn from Table 1, 
and constrained the stellar radius of each model. We assumed a metallicity ZjZQ of 0.2 for 
the SMC stars and 0.5 EMC stars, as argued in Paper 1. Following Repolust et al. (2004), we 
adopted a micro-turbulence velocity of 10 km s“^ for the models with effective temperatures 
of 36,000 K and below, and 0 km s“^ for the hotter models. 



For each star, we ran a series of grids, allowing then the values for logs', cind M to 
vary, until we had satisfactory hts (as judged by eye) between the synthetic and observed 
spectra. In a number of cases, a good £t to the strengths of the He I and He H lines 
necessitated a slight increase in the He/H ratio, as described below. The lines examined 
included the H(5, Hy, H/5, and Ho; Balmer lines, the He I A4387 singlet line, the He I A4471 
triplet line (and occasionally the He I A4713 triplet line), and the He H A4200,4542 lines. 
We examined, but did not give much weight, to the fit of the wind-sensitive He H A4686 line 
in determining the hts, but illustrate the agreement between the models and the observed 
spectra. 

As discussed both by Repolust et al. (2004) and in Paper I, FASTWIND does not 
produce a strong enough He I A4471 line for giants and supergiants of spectral type 06 and 
later. Fortunately the He I A4387 singlet line can be used for these spectral types, with 
simultaneous good fits to this line and the He H absorption. A detailed comparison of the 
output of FASTWIND and CMFGEN by Puls et al. (2005) has demonstrated that the same 
problem should be encountered with CMFGEN (since the agreement of the He I triplet lines, 
including He I A4471, in both codes is remarkable), consistent with our statement in Paper I 
that this is an example of a long-standing problem with stellar atmosphere codes (see Voels 
et al. 1989). This so-called “generalized dilution effect” is yet to be understood^. 

Below we discuss the individual stars in our sample. In determining the spectral types, 
we first consider the visual impression of the spectra, following the premises of Walborn & 
Fitzpatrick (1990). We then measured the equivalent widths of He I A4471 and He H A4542, 
the principle classification lines, and determined log IF' = logIF(4471) — logIF(4542), and 
compared that to that used to differentiate the spectral subtypes by Conti & Alschuler 
(1971). In practice, there was no disagreement between these two methods. For the earliest 
type stars, we also make reference to the new classification criteria proposed by Walborn et 
al. (2002, 2004), based upon the strengths of the N IV A4058 and N HI AA4534,42 emission 
lines and the N V A4603,19 absorption lines. Walborn et al. (2002) used this to attempt to 
break the degeneracy of the 03 spectral class and introduced types 02 (N IV>>N HI), and 
redefined 03 as “weak He I with N IV>N HI”. In Sec. 4 we attempt to provide a slightly 
more quantitative description. One needs to keep in mind, however, that there is yet to 


^Puls et al. (2005) does reveal very significant differences between the strengths of the He I singlet lines 
predicted by CMFGEN and FASTWIND in the temperature regime 36000 to 41000 (for dwarfs) and 31000- 
35000 (for supergiants). This disagreement between the behavior of He I singlets may well account for some 
of the systematics between our effective temperature scale and those of some CMFGEN studies, as well as 
the strong differences in results on similar stars found by different CMEGEN studies; see the discussion in 
Sec. 5.1. 
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be a demonstration that stars classified as “02” are necessarily hotter than stars classified 
as “03”, although this is clearly the implication of the scheme; we also examine this in 
Sec 4. Furthermore, as emphasized in Paper I, there is no a priori basis for assuming that 
the strengths of the “Of’ characteristics (He II A4686 and N III AA4634,42 emission) will 
scale with luminosity in the same way in the SMC than in the Milky Way for stars of the 
same effective temperatures, given the factor of 5 difference in metallicity between these two 
systems, and the corresponding differences in stellar winds. N III AA4634,42 emission is a 
complex NLTE effect and its size mainly due to effective temperature (Mihalas & Hummer 
1973, Taresch et ah 1997), while He H A4686 will depend upon the stellar wind properties, 
such as density and temperature. Accordingly, we will note any discrepancies between these 
traditional luminosity indicators and the actual My of the stars. 


The model fits give us a measure of the effective surface gravity g^efr which will be the 
combination of the true surface gravity ^^true offset by centrifugal acceleration due to the 
rotation of the star. This correction is estimated (in a statistical sense) as simply the square 
of the projected rotational velocity divided by the stellar radius: 


fi^true T 


usm* 


( 1 ) 


6.967? 

(Repolust et ah 2004), where the numerical factor allows the use of the usual units; i.e., 
km s“^ for usini and solar units for R. We have a good estimate of usini from comparing 
the model lines with the observed spectra, except that our spectral resolution results in 
a minimum “rotational” velocity of (generally) 110 km s“^. However, for these stars with 
minimal rotational velocities the centrifugal correction is tiny (0.01 in loggf), and insignificant 
compared to the typical fitting uncertainty of 0.1 in logg (Paper I and Sec. 5). However, for 
the fast-rotators in our sample the correction can be marginally significant, amounting to 
0.05 dex in several cases, and 0.13 dex in one case (AV 296, from Paper I). The spectroscopic 
masses should be determined from the true surface gravities; i.e., Mspect = (fi'true/fi'©)7?^- 


The final values for the fits are given in Table 5. If we were only able to determine a 
lower limit to Tgfj, then the values for logg, R, M, and the derived Mgpect are uncertain. 
Note that in Paper I the stars with only a lower limit on Rs had their stellar radii incorrectly 
listed as lower limits. 


3.1. Comments on Individual Stars 

3.1.1. SMC 


AV 111 .—This star was previously classified as 05 V by Crampton & Greasley (1982). 
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Visually, we would classify the star as an 04 V((f)) (Fig. la). This is conhrmed by our 
measurement loghF' = —0.75. The strength of He I A4471 (EW=170mA) would preclude 
it being classihed as an 03. The luminosity class “V” based upon the strong He H A4686 
absorption is consistent with the star’s absolute visual magnitude, My = —4.8. There is a 
trace of emission at N HI AA4634,42 (EW~ —250mA) leading to the “((f))” designation; 
there is no sign of N IV A4058 emission. We obtained satisfactory hts with a slightly increased 
He/H number ratio (0.15), as shown in Fig. lb. The mass-loss rate is low (< 5 x IO^^Mq 
yr“^), and only an upper limit can be established. (A value of 3 x IO^^Mq yr“^ is used in 
the htting.) The Ha prohle obtained from the ground-based data agrees well with the HST 
data; we use the former, since its S/N is greater. 

AV 435 .—This star had been previously classified as 04 V by Massey et ah (1995). Here 
we would call it an 03 V. Although He I A4471 is weakly present (Fig. 2a), its equivalent 
width is only 125mA. This is comparable to or smaller the strength seen in the original 03 V 
prototypes HD93128 and HDE 303308 (Simon et al. 1983), as well as the other star used by 
Walborn (1971a) to dehne the 03 class, HD 93250 (Kudritzki 1980). Walborn et al. (2002) 
has proposed revising the early 0-type classes based upon the strength of selective emission 
lines, e.g., N IV A4058 and N HI A4334,42. Since N IV A4058 is (weakly) present (EW=- 
80mA), with N HI emission weaker (EW> —60mA), this would also lead to an 03 V((f*)) 
classihcation according to these revised criteria. The ((f*)) simply refers to the fact that 
N IV A4058 is in emission and equal or stronger than any N HI emission. The absolute visual 
magnitude My = —5.8 would be consistent with either the “V” or “HI” luminosity class; see 
Conti (1988). The hts we obtained (Fig. 2b) are a good match to the observed spectrum. 
Again we found the agreement between the ground-based Ha prohle and that obtained with 
HST is good, and show the former in the hgure. 

AV 440 -— Previously the star had been classihed as 07 V by Garmany et al. (1987), 
but it is clear in our spectrum (Fig. 3a) that He I A4471 is somewhat stronger than He H 
A4542, and visually we would call this an 08 V star. We measure log IF' = +0.12, in accord 
with this determination. The absolute magnitude My = —4.9 is consistent with the “V” 
luminosity class suggested by the star’s He H A4686 absorption. 

In order to hnd a good ht, we needed to again slightly increase the He/H number ratio 
from the canonical 0.10 to 0.12. This allowed us to then obtain excellent matches to the He H 
lines (Fig. 3b). Although the central core of the the model’s He I A4471 line is weaker than 
observed, the EWs of this line in the adopted £t and from the star are nearly identical (0.70A 
vs. 0.78A, respectively), and we find excellent agreement for the He I A4387 line as well as 
the He I A4713 line (not shown). No ground-based Ha was obtained, and so we rely upon the 
HST observation for the mass-loss determination. Although we used a value of 1 x 10“'^ Mq 
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yr in fact we could only set an upper limit on the mass-loss rate (M < 3 x 10 yr ^). 

AV 446 -—This star was classihed as 06.5 V by Garmany et al. (1987). We measure 
log W = —0.13, consistent with this classihcation. The only spectral anomaly is that the line 
depth of He I A4471 is somewhat greater than that of He H A4542 in our spectrum (Fig. 4a), 
even though the He I line has a smaller equivalent width. This implies that the line widths 
differ. This is usually an indication of a binary (see, for example, Walborn 1973a and the 
follow-up study by Massey & Conti 1977). However, the absolute magnitude My = —4.7 is 
indicative of a single star with luminosity class “V”, consistent with the spectral appearance. 
The good agreement of the model and observed spectra (Fig. 4b) also supports a single-star 
interpretation. A comparison of the ground-based and HST Ha prohle shows reasonable 
agreement, and so we used the higher S/N ground-based data for the £t. Once again we 
could set only an upper limit on the mass-loss rate. 

AV476 .—This star was classified as an 06.5 V by Massey et al. (1995). Our higher S/N 
spectrum, shown in Fig. 5, would suggest a somewhat earlier type. We measure logW = 
—0.34, indicative of a spectral type of 05.5. He H A4686 is comparable in strength to He H 
A4542, and there is no sign of N HI AA4634,42 emission. Thus if this were a Galactic star, 
we would assign a luminosity class of “V”. The high absolute visual luminosity of this star 
My = —6.3 is more in keeping with the star being a giant. We argued in Paper I that these 
“f’ characteristics should not (and often do not) scale with luminosity the same way in the 
metal-poor SMC as they do in the Milky Way, as the He H A4686 emission that is used as 
luminosity indicator will be weaker (at the same luminosity) in the SMC as the stellar winds 
are weaker. An alternative explanation could be that this star is a binary, and the spectrum 
composite. Our fitting efforts are consistent with the latter interpretation in this case, as the 
hydrogen and He H lines have a very high radial velocity (270 km s“^) which is not shared 
by the He I velocities. There is a hint of N IV A4058 emission in the spectrum (Fig. 5), and 
we suggest that this star is an 02-3 V plus a somewhat later 0-type. We note that there is 
significant nebular contamination in the ground-based Ha profile. 


3.L2. The non-R136 LMC Stars 

Sk—67°22. —Classically, we would classify this star as 03If*/WN, where the second 
part of the nomenclature simply denotes that the emission features are comparable to what 
might be found in a Wolf-Rayet star (see Fig. 6a). Indeed, this star was one of the first 
to be described with this “intermediate” designation (Walborn 1982). Nevertheless, its 
spectrum is quite unlike that of a typical Wolf-Rayet star, in that its absorption-line spectra 
lacks P-Cygni components, and it is much more like an ordinary 0-type star than more 
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extreme “intermediate” stars, such as Br 58 discussed below. Nevertheless, the “slash” 
designation resulted in its being included in the recent compilation of Wolf-Rayet stars in 
the LMC (Breysacher et ah 1999). Were it to be classihed according to the Walborn et al. 
(2002) criteria, it would have to be considered to be an 02If*, as N IV A4058 emission is 
much stronger than the slight NIII AA4634,42 emission visible on the blue wing of He II 
A4686, and there is no hint of He I A4471 in our spectrum. However, Walborn et al. (2002) 
explicitly excluded 03If*/WN objects from consideration when discussing the 02 spectral 
class. Following in these footsteps, should we call this star an 02If*/WN? It seems to us 
that all of the so-called “03If*/WN” stars might be similarly reclassihed as 02 If*. Both 
Walborn (1982) and Conti & Bohannan (1989) have stressed the difficulty of classing such 
“intermediate” objects, and the introduction of the 02 I classihcation has added yet another 
wrinkle. Furthermore, the prototype (and until Paper I, the only member) of the 02 If* 
class is the star HD 93129A, which has relatively weak He H A4686. However, Nelan et 
al. (2004) has now split HD 93129A into two components (Am = 0.5 mag) separated by 
OOmas, which might account for the weaker emission. In Sk—67°22 the NIV A4058 emission 
is almost identical in strength to that of R136-020 (discussed in Paper I), while the He H 
A 4686 line is nearly twice as strong. Several additional examples will be discussed below, 
and a comparison of their spectra made in Sec. 4. C IV A4658 is present (EW=-300mA) in 
emission. 

Despite the strong emission, our htting was straightforward, and good hts were obtained 
(Fig. 6b). We had to increase the He/H ratio to 0.3 to obtain He lines of the right strength 
relative to H. Only a lower limit could be set on Tgfj since no He I could be detected. However, 
if the temperature were 45,000 K rather than 42,000 K we would have to signihcantly increase 
the He/H ratio even further to maintain a good £t. We note that this star has a large 
projected rotational velocity (nsini = 200 km s“^) compared to most of the other stars in 
our sample, and we offer the speculation that the chemical enrichment implied by He/H~ 0.3 
and the strong emission are both coupled to this. As Maeder & Meynet (2000) demonstrate, 
high rotational velocities can result in efficient mixing of processed material from the core of 
a star, and it can also lead to enhanced mass-loss, particularly from the poles (see Owocki 
et al. 1998). 

The radial velocity of this star is quite large, 430 km s“^, and is about 150 km“^ greater 
than the systemic velocity of the LMC (Kim et al. 1998). A similar radial velocity was seen 
in both of our blue CTIO spectrum and HSTHa spectrum; otherwise we would be inclined 
to believe this might be a binary. This star was cited by Massey (1998) as an example of an 
isolated (held) 03 star. Its high radial velocity suggests that it is a runaway star, ejected 
from its birth place. 



5'A:—65°47.—This star was classified as 04 If on the basis of the present spectrum 
(Fig. 7a) by Massey et ah (2000). Here we concur with this type. He I A4471 is only 
weakly present, with an equivalent width IV = 115mA. Given the weakness of He I A4471, 
we might be tempted to reclassify the star as 03, were it not for the fact that log W = —0.68, 
consistent with the 04 classification. Furthermore, Walborn et al. (2002) suggest that for 
an 03.5 N HI A4634,42 emission is similar in strength to N IV A4058 emission, while here 
we find N HI > N IV (although the latter is clearly present), also leading us to the 04 
subtype. The strengths of the N HI A4634,42 and He H A4686 emission lines lead to the “If’ 
luminosity class. This is consistent with the absolute magnitude My = —6.4, a value which 
is typical for an 04 If star (Conti 1988). Good fits to the lines were easily found (Fig. 7b), 
and the physical parameters are well constrained. 

LH58-496 .—This star was classified as 03-4 V by Conti et al. (1986) and Garmany et 
al. (1994). Our higher S/N data (Fig. 8a) suggests a later type, 05 V, consistent with our 
measurement of log W = —0.52. The absolute magnitude. My = —5.1, is in good agreement 
with our spectral classification of a dwarf. We obtained a good fit (Fig. 8b), except for the 
H/3 profile which is likely filled in by emission. The Ha data were obtained with EST and 
thus provides a relatively nebular-free profile. 

LH81:W28-23 .—This star was previously classified as 03 V by Massey et al. (2000) 
based on the same data as used here. Although He I A4471 is readily identified in our 
spectrum (Fig. 9a), its presence is due to the high S/N, as we measure an equivalent width 
of only 80mA. We find a value of log IF' = —1.0, well below the value —0.6 that separates 
the 04’s from the 05’s (Conti 1988). N V AA4603,19 absorption is clearly present, and N IV 
A4058 emission (EW -120mA) is less than that of N HI AA4634,42 (EW = -500mA). Thus 
under the Walborn et al. (2002) scheme, the star would be classified as “03.5 V((f+))”, with 
the “((f+))” denoting N HI emission with comparable Si IV emission. The strong He H 
A4686 absorption indicates that the star is a dwarf, consistent with its My = —5.1 value. 

Our modeling of this star was straightforward, although it quickly became clear that we 
had to significantly increase the He/H number ratio from its canonical 0.1 value to 0.25 in 
order to match the observed spectrum. The final fits are shown in Fig. 9(b). The ground- 
based Ha profile agrees well with that obtained with HST-^ accordingly, the former is shown, 
as it is considerably less noisy. 

LH90:Br58 .—Breysacher (1981) included this star in his catalog of LMC Wolf-Rayet 
stars, where it is classifed as WN5-6. The star was reclassified by Testor et al. (1993) as 
WN6-7. The lower Balmer lines are in emission (Ha) or P Cygni (Hy, H5), and He H A4200 
and A4542 are present primarily in absorption (Fig. 10), leading Massey et al. (2000) to 
argue that the star is more properly called an 03 If/WN6, although as mentioned earlier the 
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classification of such “intermediate” objects is difficult. The equivalent width of He II A4686 
is — IsA; usually —10 is taken (somewhat arbitrarily) as the dividing line between an Of star 
with strong emission line and a bona-hde Wolf-Rayet star of the WN sequence. Clearly the 
emission is considerably greater than in Sk—67°22 (Fig. 6a). Although Br 58 is a member 
of a “tight cluster” (see Fig. la of Testor et ah 1993), our 1.3” slit width and position angle 
excluded resolved neighbors. We note that there is little evidence from a visual examination 
of the spectrum that there is enhanced helium present at the surface of the star, and a more 
detailed analysis than what we can perform is needed to argue if nitrogen is enhanced or 
not. 


In any event, the strong emission defeated our efforts to obtain an acceptable £t. We did 
establish that a steep velocity law {{3 = 3.0) and high mass-loss rate (40 x yr“^) gave 

an excellent £t to the Ha profile, with a temperature of 40,000-42,000 K giving a reasonable 
match to the He H line strength if a logg = 3.5 was assumed. The high value of (3 alone 
might suggest that the star is physically more like a Wolf-Rayet star than like a normal O 
star, for which this large a value for f3 would be unexpected. Our £t is not of the same quality 
as that we achieve for other stars, as clearly the treatment of the stellar wind dominates the 
stellar parameters. If these parameters are even approximately correct, the star is quite 
luminous (Mboi = —11.1) with a radius of 30Rq and hence an inferred spectroscopic mass 
> IOOMq. An analysis that includes separate determinations of the elemental abundances 
(not possible with our current implementation of FASTWIND) is clearly warranted for this 
interesting object. We do not include this star in our analysis in Sec. 5, and therefore it has 
no effect on our derived effective temperature scale. In addition, we exclude this star from 
consideration in Sec. 4 since its physical parameters are poorly determined. 

LH90:ST2-22 .—This star was classified as 03 HI(f*) by Schild & Testor (1992), and 
reclassihed as 03 V((f)) by Massey et al. (2000) based upon the same spectral data as used 
here. The spectrum in shown in Fig.11a. He I A4471 is weakly present (EW of 140mA), with 
log W = —0.85. Should this star be instead classified as an 04? We find N IV A4058 emission 
weaker of N HI AA4634,42 (EW of —200mA versus EW=—700mA), with the presence of 
N V AA4603,19 absorption uncertain. Thus according to the criteria given by Walborn et 
al. (2002) we would classify this star as an 03.5 were we to rely upon the nitrogen emission 
lines. As to the luminosity class, the prototype 03.5 HI(f*) star Pismis 24-17 (Walborn 
et al. 2002) clearly has much stronger N V AA4603,19 absorption (compare Walborn et al. 
2002 Fig. 4 with our Fig. 11 a), although the N IV and N HI emission appears stronger 
in LH90:ST2-22 than in the 03.5 V(f-I-) star HD 93250. The appearance of He H A4686 
suggests a giant (HI) luminosity class. The absolute magnitude of LH90:ST2-22 is —6.4, 
which is more typical of an 03-4 s^pergiant (03-4 I) than a dwarf (Table 3-1 of Conti 1988). 
We adopt 03.5 Hl(f-I-) as the spectral type, where the “-I-” denotes that Si IV AA4089,4116 
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is in emission, a redundant reminder given that all stars so classified have this feature (see 
Walborn et ah 2002). 

The fitting of this star was straightforward. The widths of the Balmer line wings neces¬ 
sitated a value logg = 3.7. The fits are shown in Fig. 11b. A slightly elevated He/H number 
ratio (0.2) was needed in order to get the helium lines as strong as what was observed. 

BI 237 .—This star was classified as an 03 V by Massey et al. (1995) using older, poorer 
S/N data than what we have obtained for the present study. The spectrum is shown in 
Fig. 12a. We see that this star is somewhat earlier than that of LH90:ST2-22 (Fig. 11a) 
discussed above, with He I A4471 nearly absent. We do detect some He I A4471, but at 
a very weak level (~ 20mA), a detection possible only with our extremely high (500) S/N 
spectrum. With the EW of He H A4542 ~ 750mA we conclude logW' ~ —1.6! The nitrogen 
emission line criteria proposed by Walborn et al. (2002) would make this either an 02 V((f*)) 
[“NIV » NHI, no He I”] or an 03 V((f*)) [“NIV > or si N HI; very weak He I”], given 
that N IV A4058 emission (EW=-165mA) is stronger than the N HI AA4634,42 emission 
(EW=-90mA). 

Whether one should call the star an 02 or an 03 is of course dependent upon the in¬ 
terpretation of the non-quantitative criteria enumerated by Walborn et al. (2002). However, 
the presence of C IV A4658 emission, with an EW=-140mA, (also seen in BI 253, discussed 
below) is taken by Walborn et al. (2002) as an extreme hot 02 spectrum. We therefore 
call the star an 02 V((f*)), despite the (very) weak presence of He I A4471. The absolute 
magnitude (My = —5.4) is consistent with the dwarf designation suggested by the strong 
He H A4686 absorption feature (Fig. 12a). 

Despite the extreme spectral type, the fitting was very straightforward, with the only 
complication that the He H A4200 profile was corrupted. The fit is shown in Fig. 12b. The 
physical parameters are extremely well constrained thanks to the very weak presence of He I 
A4471. The HST Ha profile agrees well with the higher S/N ground-based data, and thus 
we used the latter. 

This star is one of the “field” early-type stars identified by Massey et al. (1995), located 
well outside the nearest known OB association (LH 88). We do note that the radial velocity of 
this star is quite high (~ 400 km s“^), compared to the systemic LMC velocity of 279 km s“^). 
Thus this field 02 star would also qualify as an OB runaway. In concert with the similar 
results for the field 03 star Sk—67°22, we believe this provides an explanation for the long¬ 
standing mystery of such early-type field stars (Massey et al. 1995; Massey 1998). 

BI 253 .—This star has many similarities to BI 237. It too was classified as 03 V by 
Massey et al. (1995), but was part of the basis for Walborn et al. (2002)’s extension of the 
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spectral sequence to 02. Based upon the spectrum we obtained as part of the present study, 
Walborn et ah (2002) propose BI 253 as the prototype for the newly dehned 02 V class. It 
is interesting to compare its spectrum (Fig. 13a) with that of BI 237 (Fig. 12a), which we 
argue above is also a member of the 02 V class. The nitrogen spectrum would suggest that 
BI 253 is somewhat hotter, as N III AA4634,42 emission is not visible in our spectrum, while 
N IV A4058 emission and N V AA4603,19 absorption are both considerably stronger than 
in BI 237. The EW of NIV A4058 is nearly three times stronger in BI 253, with the line 
having an EW of about -450mA. He I A 4471 is not detectable, although this could simply 
be due to the poorer S/N of this spectrum compared to that of BI 253 (225 vs 500). We feel 
conhdent in placing an upper limit of 25mA on its presence, requiring log W < —1.5, as the 
EW of He H A4542 is 710mA. C IV A4658 emission is present (EW=-160mA) in comparable 
intensity to that of BI 237. The star’s absolute magnitude (My = —5.5) is consistent with 
the dwarf designation suggested by the appearance of He H A4686. We describe the spectrnm 
as 02 V((f*)), consistent with Walborn et al. (2002) notation, althongh we note that the 
lack of N HI A4634,42 emission makes the use of the “f ’ designation arguable. 

The htting was against straightforward (Fig. 13b), although only a lower limit could be 
set on the effective temperature (and hence on other stellar parameters) due to the lack of 
an nnambiguous detection of He I A4471. The mass-loss rate is somewhat greater in BI 253 
than it is in BI 237. 

LH101:W3-14 -—Both LH101:W3-14 and LMC2-675 (discussed below) are stars whose 
spectra are moderately early (~ 05) according to the He I/He H ratio, bnt which show nitro¬ 
gen features typical of much earlier stars. Are these stars spectral composites, or nitrogen- 
enriched objects? We have previously shown that when a star is a spectral composite we 
can seldom match the actual He I and He H line strengths with a single model. The star 
was classified as 03 V by Testor & Niemela (1998), and then reclassihed as ON5.5 V((f)) 
by Massey et al. (2000) based on the spectrum used here. The Testor & Niemela (1998) 
spectrum of this star is clearly noisy (see their Fig. 4), and it is not clear whether or not He I 
absorption is present in their spectrnm or not. It is certainly present in ours, at a strength 
typical of an 05 star, with a value of log IF' equal to —0.49. Here we also see (Fig. 14) 
N IV A4058 and N HI AA 4634,42 emission of comparable intensity, so if we were to ignore 
the He I absorption spectrum we would conclude that the star was 03 V((f*)) according 
to the criteria of Walborn et al. (2002). N V AA4603,19 absorption is also weakly present, 
characteristic of an early spectral type. On the other hand, the presence of C HI A4650 
absorption is typical of a later spectral type (08.5+). The absolute visual magnitude of the 
star. My = —5.6, is somewhat brighter than the My = —5.2 value listed by Conti (1988) for 
an 05 dwarf, and is more consistent with a giant classihcation. The strength of He H A4686 
absorption is weaker than that of He H A4542, and so a giant designation is not preclnded. 
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Alternative, the absolute magnitude may be suggesting that the object is the composite of 
two dwarfs. 

The modeling quickly demonstrated to our satisfaction that this spectrum is composite. 
A 40,000 K model can match the line strengths of the He II lines and that of He I A4387, 
but this model’s He I A4471 line is then many times stronger than what is actually observed. 
(This goes in the opposite sense of the “dilution effect” described by Repolust et al. 2004.) 
We can achieve the correct line ratio for He I A4387 and He I A4471 by lowering the effective 
temperature to 38,000 K, but the He I line strengths are then much greater than what is 
actually observed, and the He H lines are too weak. This is consistent with a composite 
spectrum, with an 03 V dominating the nitrogen and He H spectrum, and a latter 0-type 
dominating the He I spectrum. We further note that the He I A4471 line appears to be 
broader (nsini = 130 km s“^) than the He H lines (nsini = 110 km s“^), a difference 
that is quite noticeable with our S/N and dispersion. This is consistent with the spectral 
classihcation of Walborn et al. (2002) based upon the same data. 

LH101:W3-19 .—This star was classified as 03 If by Testor & Niemela (1998) and again 
as 03 If* by Massey et al. (2000); the data for the latter are the same as what are used here. 
Using the criteria suggested by Walborn et al. (2002) the star would have to be considered 
of spectral type 02 If*, as the N IV A4058 emission (EW=-300mA) is greater than the N HI 
AA463442 emission (EW=-230mA). (Fig. 15a). The presence of weak C IV A4658 emission 
(EW=-150mA) suggests it is quite hot. N V AA4603,19 absorption and Si IV AA4089,4116 
emission are clearly visible. He H A4686 is in strong emission rather than P Cygni emission, 
attesting that it is a supergiant and not a giant. A supergiant classihcation is consistent 
with the high absolute visual magnitude. My = —7.04, found for this star. 

We have HST and ground-based data for both the blue part of the spectra and Ha. 
Nebulosity is strong around this star, and the HST data proved crucial for the Ha region. 
Although the HST data in the blue give less contaminated Balmer lines, the S/N is so much 
worse (130 vs. 420 per 1.4A) that we used the ground-based data for the htting, but checked 
that the HST data yielded similar results. 

The htting of this star was relatively straightforward, with reasonably good hts. The 
surface gravity of the star is somewhat higher than what we might naively expect for a 
supergiant, with logg = 3.9. The inferred spectroscopic mass is very high, > I 9 OM 0 , one 
of the highest known. He I A4471 may be marginally detected, with an EW of 50mA (and 
loglU' ~ ~1-1); we were unable to convince ourselves whether or not this feature was real, 
and we therefore treat the ehective temperature as a lower limit. 

LMC2-675 .—This star was classihed as “03 If-l-0” by Massey et al. (1995), who recog- 
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nized its composite nature. The spectrum of this star (Fig. 16) is very similar in appearance 
to that of LH101:W3-14 (Fig. 14), with a He I to He H ratio suggesting a spectral type 05, 
but with strong N features that would be characteristic of an earlier type. Our modeling 
conhrms that this object is a composite. No model simultaneously hts the strengths of the 
He I and He H lines. In this case we find we can match the He H lines and He I A4387 line 
with a single model of temperature 40,000 K, but that the He I A4471 line from the model 
is much too strong. Furthermore, the He H lines are shifted in radial velocity by about 
30 km s“^ relative to those of He I and H. The strength of N IV A4058 emission and N V 
A4603,19 absorption, plus the absence of any N HI A4634,42 emission would suggest that 
one of the components is an 02 giant [i.e., 02 HI(f*)]. The other component must be an 
0-type dwarf. Note that for the earliest types (02-3.5) the absolute visual magnitude is not 
a good indicator of luminosity class, as argued in Paper I, so this need not be in conflict 
with the relatively modest My = —5.0. 

3.1.3. The R136 Stars 

The R136 cluster in the LMC contains more extreme 0-type stars than are known in 
total elsewhere (Massey & Hunter 1998). For an interesting subsample of stars, the FOS 
“classihcation quality” spectra of Massey & Hunter (1998) have been supplemented by higher 
resolution STIS spectra with considerably greater S/N for the purposes of modeling. Since 
the STIS spectra cover only the wavelength region from Hy to He H A4542 (plus a separate 
exposure at Ha; see Sec. 2), we make use of the FOS spectra in determining the spectral 
type, but use only the STIS data for the modeling. The model spectra for the He H A4686 
line are, however, compared to the FOS data for this line. In presenting the blue-optical 
spectra, we have spliced in the better STIS data in the wavelength region 4310A to 4590A. 

R136-007 .—This star was classihed as 03 If*/WN6-A by Walborn & Blades (1997); 
Massey & Hunter (1998) agreed with this spectral type. The spectrum is shown in Fig. 17. 
As noted earlier, the “slash” designation simply denotes strong Of-like emission. Walborn 
et al. (2002) specihcally excluded such “slash” stars from their discussion, but to us this 
division seems arbitrary and unnecessary. The absorption spectrum is straightforwardly 
interpreted as 02-03 given the strong He H A4542 and lack of detection of He I A4471. 
Based upon the nitrogen emission spectrum, we would classify the star as 02 If*, as N IV 
A4058 emission is much stronger than N HI AA4634,42 emission. C IV A4658 emission is 
also present, suggesting this is particularly hot 02 star. The presence of N V AA4603,19 
absorption and Si IV AA4089,4116 emission is also consistent with this designation. The 
EW of He H A4686 is -7A, somewhat short of the -lOA (arbitrary) boundary between WRs 



and Of stars. The intrinsic emission will be slightly stronger, given the presence of a fainter 
companion (see below). Its absolnte visnal magnitnde My = —6.9 is similar to that of the 
02 If* star LH101:W3-19 discussed above, and is consistent with a supergiant designation. 

We had difficulty matching the depth of the hydrogen and He II lines with our modeling, 
and are forced to conclude that this star is a composite. Indeed, Massey et ah (2002) found 
light variations of this star that are indicative of eclipses. Subsequent imaging has now 
identihed a period for the system, and a radial velocity study is in progress. 

R136-014 -—This star was classihed as 04 If by Melnick (1985), under the designation 
Mk37) and as 04 If+ by Massey & Hunter (1998). Although He I A4471 is weakly present in 
our STIS spectrum (Fig. 18a, the EW is only 40mA, and it is difficult to see how it could have 
been unambiguously detected in these earlier studies. We derive a value of log IF' = —1.0, 
suggesting a type earlier than 04. Walborn et al. (2002) reclassify the star (under the name 
“MH14”) as 03.5 If*, given the equal strengths of the N IV A4058 and N HI AA4634,42 
emission lines. The He H A4686 line has an EW of -7A, and the star’s absolute magnitude. 
My = —6.4, is consistent with it being a supergiant (Conti 1988). 

The hts were relatively straightforward, and the comparison between the model and 
observations are shown in Fig. 18b. Note that although the He H A4686 was not used in the 
process, the agreement between the model and the observation is quite good. 

R136-018 .—This star was classihed as 03 HI(f*) by Massey & Hunter (1998), and this 
spectral type was repeated by Walborn et al. (2002). We show the spectrum in Fig. 19a. He I 
A4471 is present with an EW of 90mA; we measure log IF' = —0.87. NIV A4058 is present 
in emission (EW=-750mA), although NHI AA4634,42 is not; the N emission lines would 
therefore allow either an 02 or 03 classihcation according to Walborn et al. (2002), while 
the strength of He I would argue for the later type. The strength of He H A4686 absorption 
and overall weakness of the nitrogen emission spectra argues that the star is a dwarf or giant 
rather than a supergiant. Its absolute magnitude. My = —5.9, is intermediate between 
that expected for a dwarf and a supergiant (Conti 1988), and so we retain the 03 HI(f*) 
classihcation. 

The htting of this star required only two model runs to obtain an excellent ht; we show 
the agreement between the observation and model spectrum in Fig. 19b. The surface gravity 
is found to be consistent with a giant. 

R136-033 .—This star was classihed as 03 V by Massey & Hunter (1998). He I A4471 
is marginally detected, with an EW 50mA and log IF' = —1.2 (Fig. 20a). By the “old” 
criteria, we would call this spectral subtype 03. Neither N IV A4058 nor N HI AA4634,42 
emission is detectable in our FOS spectrum, and so it is difficult to know how to apply the 
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criteria suggested by Walborn et al. (2002) based upon the relative intensities of these lines. 
(We will return to issue in general in Sec. 4.) For now, we will call the star an 02-3.5 V. 
The absolute magnitude My = —5.1 is consistent with the star being a dwarf. 

The htting for this star was also straightforward, and the results are shown in Fig. 20b. 
The surface gravity for this star is found to be similar that of R136-018 (loggf = 3.75), 
despite the fact that R136-033 is a dwarf and R136-018 is a giant. We note in Sec. 5 that 
there is no clear distinction in the surface gravities of giants and dwarfs in general among 
the earlier O stars. We list its parameters in Table 5 as if the presence of He I A4471 were 
real, but possibly these should be viewed as upper limits. 


3.2. Agreement with Other Wavelength Regions 

3.2.1. Far Ultraviolet 

There are five stars in our sample for which we have FUSE data. Of these, one star is 
composite (AV 476), and we do not consider this further. For the remainder, we performed 
SEI htting of the OVI in the same manner as was carried out on the other UV data. In 
Fig. 21 we compare the hts of the OVI line to that of CIV for the star Sk—67°22. In each 
panel the red line shows the adopted Uoo=2650 km s“^, with the dotted blue lines showing 
±100 km s“^. We compare the velocities for all four stars in Table 6; clearly the agreement 
is excellent. 

We can of course go much further with both the FUSE and HST/UV data. Some 
studies of the physical properties of 0-type stars have relied solely on the ultraviolet region 
for modeling (see, for example, Martins et al. 2004, Garcia & Bianchi 2004, Bianchi & Garcia 

2002, Pauldrach et al. 2001), while other recent studies have included both the optical and 
ultraviolet regions in their modeling efforts (for instance, Crowther et al. 2002, Hillier et al. 

2003, Bouret et al. 2003). In general the UV provides the much-needed diagnostics of the 
stellar wind lines (particularly the terminal velocities) but the various lines that have been 
proposed as particularly temperature sensitive (for example, C III A1176 to C IV A1169 by 
Bouret et al. quoting S. R. Heap in preparation; Fe IV to Fe V by Hillier et al. 2003, Ar VI 
AlOOO to Ar VH A1064 by Taresch et al. 1997, etc.) are usually useful over a very limited 
temperature range, and, more importantly, good matches between the model and stellar 
spectra require determining the abundances of the particular element. The lines in the UV 
region are also not particularly sensitive to the surface gravity; see discussion in Pauldrach 
et al. (2001). We believe that there is much to be gained by fitting the UV and optical 
region, as has been shown in the past by some (see in particular Taresch et al. 1997), but 
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has generally been ignored by others. Although the present paper (and Paper I) has made 
some use of the UV (in terms of determining the terminal velocity), our study has been 
wedded to that of the optical lines, as these lines have the greatest sensitivity to loggf and 
Teff. Nevertheless, the UV region contain differing ion states of various metal lines that can 
be effectively employed. In addition, wind clumping will affect the UV and the optical in 
different ways. The most robust answers will come from fully utilizing both the optical and 
UV, and we plan to make further use of our beautiful UV spectra (both HST and FUSE) in 
a subsequent paper 


3.2.2. The NUV: He I A3187 and He H A3203 

In Paper I we (re) introduced the He I A3187 and He H A3203 lines as an interesting 
check on the models and effective temperatures we adopt. Morrison (1975) was the first to 
call attention to the usefulness of these lines, noting that the He H A3203 (n=3-5) line was 
the only accessible He H line that did not involve transitions from n = 4. The He I A3187 
line is a triplet 2^S-4^P°, similar to A4471 (2^P°-4^D). 

Of the six stars in this paper for which we have observations in the NUV region, two are 
composites. In Fig. 22 we show the spectra of the other four. Although we did not use this 
region in determining the physical parameters of these stars, in all four cases there is good 
agreement between the observed spectra and the models for He H A3203. In Paper I we also 
demonstrated good agreement for the two stars for which we had NUV data. We conclude 
that the He H A3203 line yields answers which are consistent with that of He H A4200 and 
A4542. None of the stars were sufficiently late for He I A3187 to be detected at our S/N, 
consistent with the output of the stellar atmosphere models. Since He I A4471 and, in some 
cases, even He I A4387 is measurable in all of these stars (with the possible exception of 
LH101:W3-19; see Sec. 3.1), we note that while the NUV region provides reassurance on our 
htting procedure, the additional information added is limited. 


4. The Earliest O Stars 

Although the main goal of our paper is to derive physical parameters, we have assembled 
some of the highest quality optical spectra on some of the earliest O stars known, and 
we would be remiss to not use this to comment upon the spectral classihcation scheme. 
Conti & Alschuler (1971) provided equivalent width measurements of the primary spectral 
classihcation lines He I A4471 and He H A4542 for a large sample of O stars that had been 
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classified in the “traditional way” (visually comparing the spectra of program stars to the 
spectra of spectral standards) and provided a quantitative scale which could be used to 
distinguish the spectral subtypes from one another; i.e., an 08 spectral type had a value of 
logff^' = EW (He I A4471) / EW (He H A4542) between 0.10 and 0.19, while an 08.5 star 
had a value of log W between 0.20 and 0.29. This scheme was slightly revised for the earliest 
types by Conti & Frost (1977), who used the criteria that 04 stars have log IT' < —0.60 
while in 03 stars He I A4471 is “absent”. (A summary of the loglH' appears in Table 1-3 
of Conti 1988.) Of course, this leads to the 03 classihcation as being degenerate; i.e., at a 
given surface gravity, stars with any temperature above some value might be classihed as 
“03”. Indeed, when higher signal-to-noise spectrograms were obtained by Kudritzki (1980) 
and Simon et ah (1983) of the prototype 03 stars HD 93205, HD 93128, and HD 303308 
(Walborn 1971a, 1973b), weak He I A4471 was revealed, with EWs of 75-250 mA, and 
loglH'= -1.1 to -0.5. (HD 303308, with a loglH' = —0.5, would today be classihed as an 04 
star; see Table 3 of Walborn et ah 2002. Another 03 star for which Simon et al. 1983 detected 
He I A4471 is HD 93129A, but today we recognize that this is a composite spectrum; see 
Nelan et al. 2004.) However, to date no one has actually provided a quantitative distinction 
in the He I/He H ratio between spectral types 03 and 04. 

Detection of weak He I A4471 is crucial for an accurate temperature determination using 
our method, and we can ask what is our expected detection limit? With a S/N of 400 per 
1.4A spectral resolution element, we should be able to reliably detect (at the 3a level) a line 
whose equivalent width is 3 X (1/400) X I.2A, or about O.OOOA (9mA). This is consistent 
with our measurements of the 20mA in the spectra of some stars (i.e., HI 237). For the 
R136 stars, with a S/N of 100 per 0.4A spectral resolution element, a limit of about 20mA is 
expected. In a few cases, the STIS spectra had worse S/N, and only an upper limit to He I 
A4471 can be detected, and we note such cases. What about the effects of normalization on 
such weak lines? In normalizing our spectra, we used a low order cubic spline £t through 
line-free regions. Care was taken with this procedure. Comparing the continuum levels in 
our normalized spectra in the region around He I A4471 we find that this procedure worked 
very well, with some very early stars (HI 237) the normalization proved good to 0.05%; in 
general, we the agreement was good to a few tenths of a percent. What effect will that 
have on the measured equivalent widths? If our continuum placement was in error by 0.5% 
(which, we believe, would be an extreme case), then we would measure a 20mA line as having 
an EW of 28mA. Thus we expect that any error introduced by our normalization will have 
an effect comparable to that of the large but finite S/N of the spectra. 

Walborn et al. (2002) attempted to reduce the 03 degeneracy problem by using the N IV 
A 4058 and N HI A4634,42 emission lines to extend the spectral classification to spectral type 
02 (“N IV >>> N HI”) with the introduction of an intermediate type (03.5) for which N IV 
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~ N III. (The 02 class would then be the degenerate subtype.) The classihcation criteria 
are not quantitative, and the system still relies upon statements such as “very weak He I” 
or “no He I” as secondary criteria. The criticism has been made (both here and in Paper I) 
that relying upon the relative strengths of the optical nitrogen lines lacks a solid theoretical 
underpinning; i.e., although a unique spectral subclass may be dehned, it is not clear that 
a star’s effective temperature is the primary thing that distinguishes 02’s from OS’s in this 
scheme. The mechanism for N HI emission in O stars has been well described by Mihalas 
& Hummer (1973), who demonstrated the N HI A4634,42 lines will come into emission even 
in a static, non-extended atmosphere if the effective temperature is sufficiently high, due 
to a complicated NLTE effect. Thus the presence of N HI A4634,42 in dwarfs and giants 
should be related to effective temperature. However, for stars with signihcant mass-loss rates, 
an additional process known as the Swings mechanism (Swings 1948) will come into play, 
enhancing the N HI emission (see discussion in Mihalas & Hummer 1973). No such detailed 
analysis exists for the N IV A4058 line, although Taresch et ah (1997) show its behavior as 
a function of effective temperature holding the surface gravity constant and maintaining the 
same Ha prohle by slightly varying the mass-loss rate. Indeed, the Taresch et al. (1997) 
analysis of HD 93129A showed the power of using the N HI and N IV emission lines as 
a constraint on the effective temperature, i/other parameters (such as logg and M) were 
constrained by other observations. 

In order to provide a more quantitative assessment of the classihcation criteria, in Ta¬ 
ble 7 we provide for the hrst time equivalent width measurements and ratios for the He I 
A4471/He H A4542 absorption lines, and N HI A4634,42 / N IV A4058 emission lines, for 
a sample of the earliest 0-type stars. We have classihed the stars using the Walborn et al. 
(2002) criteria; otherwise, all of the stars earlier than 04 would simply be described as “03”. 
Some of the dwarfs in Paper I had no detected N HI or N IV emission and were simply called 
“03 V”; here we “revised” the type to “02-3.5 V”. (The need to detect these lines, which 
are weak or non-existent in dwarfs, particularly at low metallicities, is an obvious drawback 
to the new classihcation scheme.) 

First, let us note that the dividing line between 03’s and 04’s occurs at roughly a 
logW' = —0.8, and we would propose this as an improvement over the old description 
that He I be “absent”. However, even a casual inspection of Table 7 shows the ehect that 
mass-loss has on the derived ehective temperature for a given He I/He H ratio. LH101:W3- 
24, an 03 V((f-|-)) star, and HI 237, an 02 V((f*)) star, both have similar Tes and loggf, 
although their He I/He H ratios diher by a factor of 6! This can be readily attributed to the 
signihcantly greater mass-loss rate of HI 237 due to its larger radius, mass, and luminosity. 
(Recall from the discussion in Paper I and here that the ehect of stellar winds will be to 
decrease the derived ehective temperature for a given Hel/Hell ratio, in part because of the 
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filling in of Hel due to emission produced in the wind.) 

In Fig. 23a, we show the relationship between the Hel/Hell and NIII/NIV ratios. There 
is a reasonable correlation in this hgure, in that stars of a given luminosity class which have 
a small Hel/Hell ratio (right side of the diagram) also tend to show a smaller NIII/NIV 
emission line ratio (top of the diagram). The correlation is quite good for the dwarfs (hlled 
circles), with giants (triangles) and supergiants (open circles) exhibiting more scatter. To us 
this is as expected. Although Hel absorption is affected by stellar winds (by the hlling in of 
emission), the strength of the nitrogen emission lines are likely to be more affected by mass 
loss (due to, for instance, the Swings mechanism, as discussed above). 

We address this in a more quantitative way in Fig. 24 where we show the correlation 
of effective temperature with (a) the Hel/Hell ratio and (b) the NIII/NIV ratios. Again, 
this considers only the earliest-type O stars. For the dwarfs (hlled circles) there are very 
large changes in both Hel/Hell and NIII/NIV with little change in effective temperature. 
In other words, the extension of the spectral classihcation through 02-03.5 for the dwarfs 
is unwarranted in terms of effective temperature. For the giants (denoted by triangles) 
and supergiants (open circles) there is more of a correlation of effective temperature with 
either spectroscopic criteria, but neither is very good. Consider, for example, the star with 
the hottest temperature in our sample, LH64-16. An independent analysis (using the same 
optical data, however) yielded a similar effective temperature (Walborn et al. 2004) as what 
we derived in Paper I. However, the giant R136-018 has a NHI to NIV emission ratio which 
is at least as small, but has nearly the coolest effective temperature for any of the giants with 
a measured NIII/NIV ratio in our sample! (Recall from the discussion above that the lack 
of NHI emission in the spectrum of R136-018 would allow either an 02 or 03 designation; 
the argument for the 03 classihcation came about from the strength of He I absorption.) 
Of course, we argue in Paper I that the star LH64-16 is likely the result of binary evolution, 
and should not be considered a prototype of a new type (cf. Walborn et al. 2002). 

The degeneracy of Teff with log W is somewhat surprising, since afterall this is the basis 
for the original MK types (although not the extension to 02-3.5 described by Walborn et 
al. 2000). By examining the large grid of models computed as part of the Puls et al. (2005) 
study, we came to understand that the effect is due primarily to the dependence of log W 
on surface gravity for stars of the hottest temperatures. This is demonstrated in Fig. 25(a). 
There is also a signihcant dependence of log W upon the mass-loss rate for the hottest stars 
with the highest rates, and highest temperatures (Fig. 25b). 

It is also worth noting that the spectroscopic luminosity criteria proposed for the earliest 
types (based upon the strength of N HI AA4634,42 and He H A4686) do not prove to be a 
reliable indicator of the effective surface gravities. Although the supergiants in Table 7 do 



25 


have lower surface gravities than the dwarfs, one cannot differentiate the dwarfs and giants 
on the basis of logg. We return to this point in Sec. 5. 

Thus in terms of the physical properties of early-type O stars, and in particular the 
effective temperatures, we are forced to conclude that there is little benefit to the extension 
to spectral type 02 proposed by Walborn et ah (2002). For a star of the same effective 
temperature and similar surface gravity either the Hel/Hell ratio or the NIII/NIV ratio can 
vary by essentially the full range found between 02 and 03.5 (for example, BI 237 and 
LH101:W3-24). We believe that careful modeling of the N III and N IV emission lines, in 
concert with other lines will provide useful diagnostics (following Taresch et ah 1997), and 
modifications to FASTWIND to allow this are planned. However, our work here suggests 
that it is naive to expect that any one line ratio (N III/N IV, or even He I/He H) gives insight 
into the effective temperature of the star. Although the earliest 0-type stars (what would have 
classically been called “03”) contain stars with a significant range of effective temperatures, 
no one spectroscopic diagnostic (such as the He I/He H or N HI/NIV line ratios) provides a 
good clue as to the effective temperature without knowledge of the other physical parameters 
(mass-loss rate, surface gravity). This general principle has been previously emphasized by 
others, notably Sellmaier et ah (1993). 


5. Results 

With the work described here and Paper I, we have attempted to model the spectra 
of 40 0-type stars in the Magellanic Clouds. We succeeded for 13 SMC stars (8 dwarfs, 2 
giants, and 3 supergiants) and 20 LMC stars (9 dwarfs, 4 giants, and 7 supergiants), with 
the other 7 stars showing signs of composite spectra^. This sample is now large enough for us 
to determine preliminary effective temperature scales for O stars in the SMC and LMC, and 
compare these to that of the Milky Way, using the sample of stars analyzed by by Repolust 
et ah (2004) using the same methods and models atmospheres. We also wish to make good 
on our promise in Paper I to examine the Wind-Momentum Luminosity relation from these 
new data, and, hnally, to compare the spectroscopically derived stellar masses with those 
derived from stellar evolutionary tracks. 

In Table 9 we summarize the derived parameters for our hnal sample of 33 stars. In 
what follows we assume the htting errors as quoted in Paper I, namely an uncertainty of 


"‘For one of these stars, AV469, Evans et al. (2004), obtained a fit that they deemed satisfactory. Possibly 
we were too fussy in being dissatisified with our fit, and declaring this a composite, but a radial velocity 
study is underway. 
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about lOOOK (2-3%) in Tes, 0.1 dex in logg, and 20% in M. (The later assumes that /? is 
the same for all stars and is precisely known; the actual error on M may be signihcantly 
higher, depending upon the validity of this assumption.®) In general our values for v^o are 
good to 100 km s“^ (5%). Using the propagation of errors analysis by Repolust et al. (2004) 
[see their equation 8] we would then expect the uncertainty in the derived stellar radius R to 
correspond to A logi? ~ 0.03, or about 7%, where we have allowed for a 0.1 mag uncertainty 
in My. (The uncertainty in the radii of Galactic stars is about twice as great, given the 
much greater uncertainty in the distances and hence a larger uncertainty in My.) The total 
luminosity of the star is uncertain by about 0.12 mag in Mboi, or 0.05 in logL/L©. 

For the stars with only lower limits on we expect that the values for the stellar 
radii, Umom, and Mspect are all approximately correct; nevertheless, we do not use them in 
the analysis, in order not to mix good values with bad. 

In Fig. 26 we show the location of our stars in the H-R diagram, where we have over¬ 
plotted the evolutionary tracks of the Geneva group (i.e., Gharbonnel et al. 1993 for the 
SMG, and Schaerer et al. 1993 for the LMG). For simplicity we show only the H-burning 
part of the tracks. Since the time of these calculations various improvements have been made 
in the evolutionary code, such as the inclusion of the Vink et al. (2000, 2001) prescription for 
mass-loss, improved opacities, and, most importantly, the inclusion of the effects of rotation. 
However, a full set of these tracks are not available for low metallicities. We do show the 
new tracks (for an initial rotation velocity of 300 km s“^) for initial masses of OOMq and 
4 OM 0 from Meynet & Maeder (2004); we are grateful to Georges Meynet for making these 
tracks available to us. 


5.1. Effective Temperatures, Surface Gravities, and Bolometric Corrections 

In Fig. 27 we show the effective temperatures as a function of spectral subtype for the 
complete data set: SMG stars are shown in green, the LMG in red, and the Milky Way 
in black, with different symbols representing the various luminosity classes. Due to the 
constraints of limited observing time, and our emphasis on trying to understand the effect 
that metallicity has on the physical parameters of stars, onr LMG sample (intermediate 


5A range of 0.7 to 1.3 is reasonable for 0-type stars. Our inability to derive values for /? when Ha 
is in absorption produce a typical error of ±0.1 dex (25%), with a maximum error of ±0.2 dex, which is 
considerable (65%). We can safely say that at worse the values of M are uncertain by a factor of two if Ha 
is in absorption; however a 20-30% error is more typical of what we expect. This error is comparable to the 
uncertainity we derive in the fits by varying M while holding (3 constant. 



metallicity) is incomplete, but biased towards the earliest spectral types, while the SMC 
stars (low metallicity) was chosen to cover the full range of spectral subtypes for direct 
comparison with the Milky Way sample of Repolust et ah (2004). 

First, we can see from Fig. 27 that both in the Milky Way and in the SMC that the 
effective temperatures of supergiants are 3000 to 4000 K cooler than dwarfs of the same 
spectral types for early O stars (04-06). For stars of spectral type 08 and later, there is 
essentially no difference in the effective temperatures of dwarfs and supergiants. In both 
galaxies the giants (denoted by triangles) appear to follow the same sequence as the dwarfs 
(hlled circles). For the early O stars (04-06) SMC stars are roughly 4000 K hotter than 
their Galactic counterparts. By spectral type 08 the difference is about 2000 K, and by 09.5 
the data are consistent with no difference. Thus, by coincidence the effective temperature 
scale of SMC supergiants is very similar to that of Milky Way dwarfs. 

This is consistent with what we expect by way of the effects that stellar winds will 
have on the effective temperature scale of stars: at higher temperatures, stars with high 
luminosities (supergiants) will have higher mass-loss rates than stars of lower luminosity 
(dwarfs). At higher metallicities (Milky Way) the mass-loss rates will be higher than at 
lower metallicities (SMC). By the later spectral types (O8-B0) the differences become minor. 
While this is expected for dwarfs, as these stars have progressively weaker stellar winds with 
spectral types (as the luminosity decreases with later spectral type as well), it is somewhat 
surprising for the supergiants, as even BO-2 I stars can still have very high mass-loss rates. 
We can offer a few possibilities. Perhaps the effect is caused by the fact that the continuum 
edges of H and He I are becoming increasingly important at cooler temperatures. These 
edges would do their own “blanketing”, and block out signihcant amounts of the flux. Also 
at cooler temperatures the peak flux occurs at longer wavelengths, which we expect to then 
result in a decrease in the effects of wind blanketing. But, in all this we should keep in 
mind that we have included only one SMC late 0-type supergiant in our analysis, and that 
possibly additional data will invalidate this effect. 

The LMC data on the earliest spectral types (02-03.5) emphasizes again how diverse 
a group of stars these are. Nevertheless, some clear patters emerge, and in particular the 
fact that the supergiants are signihcantly (~ 5000 K) cooler than dwarfs of similar spectral 
types. 

With these trends in mind, we provide in Table 8 a provisional effective temperature 
scale for the three galaxies. Owing to the lack of data points for intermediate and late LMC 
O stars, we have adopted values for the LMC which are intermediate for that of the SMC 
and Milky Way; what limited data we do have (i.e., spectral type 05 in Fig. 27) supports 
this. We emphasize that this scale is not the hnal word on the subject; indeed, we were in 



a quandary as to what effective temperatures to assign for the 03-04 class where the SMC 
dwarfs in our sample have a lower effective temperature scale than the 05 Vs; we chose to 
assign a higher temperature than our analysis of AV 177 and AV 378 would indicate. The 
effective temperature scales can be made more trustworthy by observations and analysis of 
LMC stars of all luminosity classes at intermediate and late 0-type, and SMC stars of the 
earliest types, as well as SMC late-0 supergiants. Another priority but much needed would 
be studies of stars of additional O and early-B type in the Milky Way to complement the 
efforts of Repolust et ah (2004). 

We present the adopted temperature scale in Fig. 28(a). We include for comparison the 
Vacca et al. (1996) scale for Galactic stars, which is much hotter; the scales are shown in 
comparison to the data in Figs. 28(b) and (c). 

We were initially surprised to hnd that the effective temperature scales for giants and 
dwarfs were indistinguishable, but an inspection of the surface gravities in Table 9 reveals 
the reason; in general, there is no difference in the surface gravities derived for these stars, 
except for the two latest types. In our sample this might be due to the fact that all but 
these two giants are of early type (02-3.5) where, as we have previously emphasized, the 
physical properties are not obviously correlated with the details of the spectral properties. 
However, even for the Galactic stars studied by Repolust et al. (2004) there is only a modest 
difference found in the surface gravities of dwarfs and giants, with averages of 3.7 (dwarfs) 
and 3.6 (giants) found from their Table 1. Conti & Alschuler (1971) note that in the original 
MK system luminosity classes were dehned for O stars only for the 09 subclass, and this 
was based upon the strength of ratio of Si IV A4089,19 to He I A4143 (see Morgan et al. 
1943). Subsequently, Walborn (1971b) used the N HI A4634,42 emission lines and the He H 
A4686 line (either absorption or emission) to establish luminosity classihcation criteria. As 
we have previously noted, these lines are not gravity-sensitive per se but rather may reflect 
the effects of surface abundances, effective temperatures, and what is happening in the stellar 
winds (both mass-loss rate and density/temperature prohles) in a complicated manner. It 
is clear that the supergiants (both in the Milky Way and in the Magellanic Clouds) have 
lower surface gravities than do the stars spectroscopically classihed as dwarfs, but there 
is so far no evidence to support that the early- or mid-type giants (dehned primarily as 
having “weakened ” He H A4686 absorption) as having lower surface gravity than the dwarfs. 
Theoretical modeling that includes the N HI A4634,42 line will help address this important 
issue. 

However, for now let us note a particularly egregious example, namely the star LH64-16. 
Walborn et al. (2002, 2004) make a point of declaring this the prototype of the 02 giant 
class. Yet, compare its physical properties to the 02 dwarfs HI 253 and HI 237. LH64-16 may 



be hotter, but the three stars have the same surface gravities, similar stellar radii, mass-loss 
rates, and bolometric luminosities. The spectral appearance of LH64-16 differs from that 
of BI 253 and BI 237 primarily in the fact that the He II A4686 shows more of a P Cygni 
prohle in LH64-16, plus the nitrogen lines (N IV and N III emission, and N V absorption) are 
signihcantly stronger than in BI 237 and BI 253. The only physical difference we can discern 
between the two stars is one of chemical abundances (Paper I and Walborn et al. 2004). 
In Paper I we speculated that perhaps this star is the product of binary evolution. In any 
event, the luminosity criteria used by Walborn et al. (2002) in this case do not seem to be 
tied into the properties one usually uses to distinguish amongst the second MK dimension, 
such as actual luminosity. 

How does our effective temperature scale compare to that from other recent studies 
using line-blanketed models? We give the comparison in Fig. 29. First, let us consider the 
supergiants (Fig. 29a). The overall agreement is fairly modest. In particular, the results 
of WM-BASIC modeling of the UV spectra of Galactic supergiants by Bianchi & Garcia 
(2002) and Garcia & Bianchi (2004) [shown by black squares] give effective temperatures 
considerably lower than those indicated by the Repolust et al. (2003) data on which our 
scale is based. The FASTWIND study of Gyg OB2 stars by Herrero et al. (2002) [black open 
circles] is in better agreement, unsurprising since the same models and methodology was 
used. Even so, their data would suggest our Galactic scale should be hotter for the earliest 
supergiants and cooler for the latest than our scale. As for the Magellanic Gloud data, the 
results are quite mixed, and there is not good internal agreement of the various GMFGEN 
studies. The Growther et al. (2002) and Hillier et al. (2003) studies of several LMG and SMG 
supergiants would indicated lower temperatures than what we find (in general), while the 
Evans et al. (2004) analysis of an SMG 08.5 supergiant suggests a somewhat higher effective 
temperature than that of the 07 supergiants in the other studies. As discussed earlier. Puls 
et al. (2005) have found signihcant difference for the He I singlets produced by FASTWIND 
and GMFGEN in certain temperature regimes. Until this matter is resolved and the same 
data re-analyzed, it is hard to know what to make of these differences. 

For the dwarfs and giants (Fig. 29b) we again see that the FASTWIND modeling indi¬ 
cates higher effective temperatures than that obtained by the WM-BASIG modeling of the 
UV spectra by Bianchi & Garcia (2002) and Garcia & Bianchi (2004), who hnd temperatures 
that are about 4000 K cooler than the mean relationship we derive from the Repolust et al. 
(2004) data. (The single giant modeled with FASTWIND by Herrero et al. 2000 agrees with 
this mean Galactic relationship.) The WM-BASIG modeling by Bianchi & Garcia (2002) 
and Garcia & Bianchi (2004) involves careful fitting of the absolute strengths of such lines as 
G IV A1169, G HI A1176, P V AA1118,1128, Si IV AA1123,1128, etc. In Paper I we note that 
the flux distributions of FASTWIND and WM-BASIG models agree well despite the former’s 
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use of approximate line-blanketing and blocking; this is now demonstrated at length by Puls 
et al. (2005). We do not have a ready explanation for the systematic difference apparent be¬ 
tween the parameters derived by the WM-BASIC and FASTWIND modeling. However, we 
believe that the hrst step in investigating this is to use FASTWIND modeling on the optical 
spectra of the same stars studied by Bianchi & Garcia (2002) and Garcia & Bianchi (2004), 
and, similarly, to use WM-BASIG to model the UV spectra of the same stars studied by 
Repolust et al. (2004) using WM-BASIG and in order to understand the specihc differences. 
We have earlier alluded to the need for such a complete modeling effort. 

In our comparison of the dwarfs and giants we inclnde the results of the GMFGEN 
modeling of Martins et al. (2002), which, for the Galactic stars, indicate a hi^'/ier effective 
temperatnre scale than our own, in contrast to the WM-BASIG resnlts. The Martins et al. 
(2002) study is theoretical in the sense that no actnal stars were used; instead, the synthetic 
He I A4471 to He H A4542 eqnivalent width ratios were nsed to assign “spectral types” to 
the models, which were rnn with reasonable mass-loss rates. We note that this procednre 
avoids the potential modeling problem involving the He I singlets (Pnls et al. 2005) discussed 
earlier. 

For the SMG dwarfs modeled by Martins et al. (2004) with GMFGEN there is very 
good agreement with onr scale, despite the fact that their observations were all obtained in 
the UV. Of course, this again means that the GMFGEN modeling was nnaffected by any 
potential problem with the He I singlets. We do note in passing an intrinsic nncertainty of 
the placement of these stars in our diagram; Martins et al. (2004) estimated the spectral 
types based npon the optical spectra prodnced by the models, since they had no optical 
observations. 

However, there are other stndies of SMG dwarfs and giants which are at variance with 
onr results. First, there is AV 69, an OG7.5 HI((f)) star whose UV and optical spectra were 
analyzed by Hillier et al. (2003) [denoted by the green star in Fig. 29b], End a low effective 
temperatnre compared to onr scale. More distnrbing, at first blush, is the discrepancy 
between our effective temperature scale for SMG dwarfs and the analysis of stars in the 
NGG 346 H H region using TLUSTY and GMFGEN modeling with apparently similar results 
from the two codes (Bouret et al. 2003)®. However, we will note that NGG 346 is the 
strongest Ha sonrce in the SMG, and that the data used by Bouret et al. (2003) (and Heap 
et al. 2005) was obtained with short-slit echelle data, which limited their ability to do sky 


®In a recent preprint Heap et al. (2005) use the same optical data and derive identical answers from their 
modeling. The resulting effective temperature scale is considerably lower than ours (6000 K) for the early 
SMC O stars. 



31 


(nebular) subtraction^. 

However, since the Bournet et al. (2003) study does involve stars in a single cluster, it 
does give us a chance to allow stellar evolutionary theory to weigh in on the issue of the 
effective temperature scales. Using what became the Chlebowski & Garmany (1991) effective 
temperature scale, Massey et al. (1989) derived an H-R diagram for this cluster which was 
highly coeval, with an age of 2-4 Myr and a typical age spread of 1 Myr (see also Massey 
1998, 2003). The H-R diagram shown by Bouret et al. (2003) [their Fig. 12] implies an age 
spread of > 7 Myr rather than 1 Myr®. Even excluding the oldest star (NGC346-MP12, 
a nitrogen rich 09.5-BO V star which Walborn et al. 2000 suggest is not a member of the 
cluster), the age range is 5 Myr. Furthermore, although Bouret et al. (2003) fail to comment 
on this point, the 02 HI star would have the youngest age (0.5 Myr); its location in the H-R 
diagram would be consistent with it being a dwarf and not a giant. 

We compare the H-R diagram derived from their values to one which we construct 
using our new effective temperature calibration in Fig. 30. For consistency we have used the 
reddenings and spectral types adopted by Bouret et al. (2003) in making this comparison®. 
The effective temperature scale we adopt here lead to more consistent ages for the stars in 
the cluster. While this doesn’t prove our scale is right, and the Bouret et al. (2003) values 
wrong, it does emphasize an often overlooked point, namely that changes in the effective 


^The version of the data shown by Walborn et al. (2000) lacked nebular subtraction, and made no 
correction for the moonlight continuum contamination (although the solar line spectrum was removed in 
an ad hoc manner). This was partially corrected by a re-reduction of the data for the versions used by 
Hillier et al. (2003), Bouret et al. (2003), and Heap et al. (2005), although no correction could be made for 
NGC 346-368 or NGC346-487 since there weren’t enough sky pixels. We are indepted to Ghris Evans for 
correspondance on this subject. We have recently obtained our own high S/N optical spectra of the NGG 
346 stars, and an analysis is underway. 

®Note that although Bouret et al. (2003) incorrectly attribute the source of the photometry used in 
constructing their H-R diagram to Walborn et al. (2000), it was actually that of Massey et al. (1989); i.e., 
both studies started with the same UBV values, and thus the differences in the H-R diagrams due purely to 
the different effective temperature scales and slight differences in the treatment of reddenings. See below for 
more on the latter point. 

®It should be noted that the intrinsic colors used by Bouret et al. (2003) lead to low values for the 
reddenings of the earliest stars, E{B — V) = 0.05 to 0.12. Massey et al. (1989) find a somewhat greater color 
excess (0.09 to 0.15). By comparison to other early-type stars in the SMC, the reddening of the early O stars 
in NGC 346 is average or slightly high, as might be expected; Massey et al. (1995) find an average reddening 
of 0.09 for the SMC. Were the Bouret et al. (2003) intrinsic colors correct, then many of the early-type stars 
in the SMC would have reddenings less than the expected foreground reddening to the SMC (Schwering & 
Israel 1991, Larsen et al. 2000). The Bouret et al. (2003) values also lead to a progression in reddenings with 
spectral type. 
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temperature scale of O stars do have implications in the interpretation of H-R diagrams and 
star-formation in clusters; see discussion in Hanson (2003) and Massey (2003). We plan a 
re-analysis of the NGC 346 stars using FASTWIND using high S/N spectra with good sky 
(nebular) subtraction. 

In Fig. 31 we show the BCs as a function of effective temperatures. As was the case 
with unblanketed models, there is no apparent difference with surface gravity, as is shown 
by the small scatter. There is a slight shift with metallicities, with stars of low metallicity 
(SMC) having a BC which is perhaps 3-4% more negative (i.e., more signihcant) than the 
higher metallicity (Milky Way) stars. We hnd a relationship 

BC = -6.90 X logTeff + 27.99 (2) 

£t the data with an RMS of 0.03 mag. This is very similar to that fonnd by Vacca et al. 
(1996) considering unblanketed models. 


5.2. The Wind Momentum-Luminosity Relationship 

Kndritzki et al. (1995) introduced the Wind-Momentum Luminosity Relation (WLR), 
where 

oc (3) 

with X = l/oeff = l/(« — <^), (Puls et al. 1996). The force multipliers a and S have been 
described by Kudritzki et al. (1989): a would equal 1 in the case that only optically thick 
lines contributed to the line acceleration force, and would equal 0 if only optically thin lines 
contributed. Typical values are 0.5 to 0.7 (Kudritzki 2002). The parameter 6 describes the 
ionization balance of the wind, and typically has a value between 0.0 and 0.2 (Kudritzki 
2002 ). Since the value of is expected to have only a weak dependence upon the effective 
temperature and metallicity (and in a way that can be theoretically predicted from radiative- 
driven wind theory; see Puls et al. 2000, Kudritzki 2002, as well as Vink et al. 2000, 2001, 
but see also Martins et al. 2004 for counter-examples) one can in fact use Eq. 3 to find 
the distances to galaxies using basic physics combined with quantitative spectroscopy of the 
bright supergiants. 

In Fig. 32(a) we show the WLR for the stars in our sample in comparison to the Galactic 
sample studied by Repolust et al. (2004). Although there is considerable scatter, it is clear 
that the SMC stars show a lower value for Umom than do the Galactic stars, with the LMC 
stars being somewhat intermediate. A more subtle effect is that the supergiants of a given 
galaxy seem to lie somewhat higher than do the dwarfs and giants. 
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We can reduce the scatter in this diagram dramatically by correcting the mass-loss rates 
given in Table 9 for the “dumpiness” of the stellar wind. Puls et ah (2003) and Repolust et 
ah (2004) argue that any kind of instability requires some time to grow, and thus the lower 
part of the wind should be minimally affected by clumping. Thus, the only stars whose Ha 
prohle should be affected are those which for which the profiles are formed in a major volume 
of the wind, i.e., the stars with Ha in emission. If Ha is in absorption, this immediately 
suggests that the only wind contribution is from layers very close to the sonic point, and thus 
absorption type prohles should be little affected by wind clumping. If radiatively-driven wind 
theory is correct, then the WLR should be independent of luminosity class, and Repolust 
et ah (2004) derived a numerical correction factor of 0.44 by forcing Galactic supergiants to 
follow the same WLR as do the dwarfs and giants. We have applied the same correction 
factor to our data for the stars with Ha in emission; only three LMC supergiants are affected, 
if we discount the stars for which only upper limits have been determined. It is not clear 
what metallicity dependence, if any, there might be in this correction factor, but since only 
the LMC data are affected (other than upper limits), this is probably safe. We show the 
revised plot in Fig. 32(b). Indeed the supergiants in the LMC are in better accord with the 
dwarfs, just as Repolust et ah (2004) found for the Galactic stars. 

In general, the SMC stars have a considerably smaller i^mom than do the Galactic stars. 
What would we expect on theoretical grounds? Vink et ah (2001) calculates that for 0-type 
stars, M oc if the Leitherer et al. (1992) relation v^o oc is taken into account. We 

can therefore expect that 

D^^Z) (X ZO- 82 . 

If a remained constant, we would then expect that the slopes in Fig. 32(b) would be the 
same for all three galaxies, but that the intercepts would be -0.25 dex lower for the LMC 
than for the Milky Way, and -0.57 dex lower for the SMC than for the Milky Way. We 
show these relationships in Fig. 32(b), where we have adopted the Vink et al. (2001) slope 
and intercept as given by Repolust et al. (2004) for the Milky Way stars. We consider the 
agreement excellent in this diagram: the plot shows a clear effect between the Galaxy and 
the LMC along the lines prediced by theory, and if the upper limits for the SMC stars are 
close to the true values, then the SMC also seems to agree well, with the exception of two 
objects. The analysis of the UV spectra we have planned should provide better constrains 
in the cases where Ha provides upper limits only. Kudritzki (2002) has performed careful 
calculations of the effects of metallicity and effective temperature on the force multipliers, 
and his calculations show that at an SMC-like metallicity that a will be 3% (Teff ~ 40000 
K) and 15% (Teg ~ 50000 K) larger at an SMC-like metallicity than in the Milky Way. We 
would thus expect x to be somewhat smaller for the SMC, and the relationship slightly more 
shallow, which is certainly not excluded by our data. 
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In evaluating Fig. 32 one should keep in mind the typical errors discussed above, i.e., 
an uncertainty of 0.05 dex in \ogL/L q for the Magellanic Cloud objects, and of 0.15 for the 
Galactic objects. The error in logDmom is about 0.15 dex for both the Magellanic Cloud 
and Galactic objects. We show the typical error bars in the hgure. Most of the stars follow 
the relationship fine; the notable exception is the star AV 14, which shows an upper limit of 
-Dmom = 28.15 at logL/L 0 = 5.85. 

5.3. Comparison of Spectroscopic and Evolutionary Masses 

The analyses of our stars have yielded values for the “spectroscopic mass”, M^pect = 
{.Qtrue!9®)^ given in Table 9. We remind the reader that these ^^true values have resulted 
from a modest correction for centrifugal acceleration to the measured values obtained 
from the model hts. It is of interest to compare these values with the mass Mevoi which we 
derive from stellar evolutionary models based upon a star’s logL and Tefj. This comparison 
is shown in Fig. 33. We have included error bars, with the uncertainty in the evolutionary 
mass assumed to be due purely to the uncertainty in logL/LQ (i.e., A log Mevoi = —0.2A 
Mboi; see equation 4 of Massey 1998). These evolutionary masses were derived from the 
“standard” (non-rotating) models using the older opacities; we return to this point below. 

We see that most of the stars in our sample cluster around the mean relationship 
Mspect ~ Mevoi, but that both samples contain a fair number of objects in which the evo¬ 
lutionary mass which is considerably greater than the spectroscopic mass. This mass dis¬ 
crepancy was hrst found by Groenewegen et al. (1989) and studied extensively by Herrero et 
ah (1992) for Galactic stars. This problem was investigated by many authors; comparison 
with the masses determined from binaries have tended to support the evolutionary masses 
(Burkholder et al. 1997; Massey et al. 2002). The lowering of the effective temperatures for 
Galactic supergiants has resulted in decreasing, if not eliminating, the discrepancy for Galac¬ 
tic stars (Herrero 2003, Repolust et al. 2004). The expected reason is two-fold: hrst, lower 
effective temperatures imply a smaller luminosity of a star, and hence the deduced evolution¬ 
ary mass will be less. Secondly, the new line-blanketed models resulted in larger photospheric 
radiation pressure, so a higher surface gravity is needed to reproduce the Stark-broadened 
wings of the Balmer lines. 

Although the Galactic mass discrepancy was considerably worse for Galactic supergiants 
than dwarfs, we see that the discrepant stars in Fig. 33 include both dwarfs and giants, and 
an inspection of Table 9 suggests that the problem stars all have Teg > 45000K, and indeed 
all stars this hot show a signihcant mass discrepancy. We indicate these stars with hlled 
symbols in the hgure. (We note that the two stars with Tefj > 45000K in the Galactic 
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sample of Repolust et al. 2004 did not show a similar problem.) 

In determining the evolutionary masses we have relied upon the older tracks of Char- 
bonnel et al. (1993) for the SMC and Schaerer et al. (1993) for the LMC. The effects of newer 
tracks, including the effects of an initial rotation speed of 300 km s“^, is shown in Fig. 26 by 
the dotted lines for the GOMq and 40Mq tracks. We see that the inclusion of rotation will 
have a modest effect on the deduced masses compared to the non-rotating models. For the 
hotter stars the effect will be to make the evolutionary masses lower; at cooler temperatures 
(on the main-sequence) the effect is in the other direction. 

For some of the discrepant stars in Fig. 33 use of the newer evolutionary models could po¬ 
tentially bring the evolutionary masses into closer agreement with the spectroscopic masses. 
However, the effect is small compared to the high temperature mass discrepancy we note: 
for stars near the ZAMS, the difference is about 0.25 mag, which is equivalent to 0.05 dex 
in logM, or about 10% 

It is possible that the objects with large mass discrepancies in Fig. 33 represent the re¬ 
sults of binary evolution. A good candidate is LH64-16, discussed in Paper I. This “ON2 III” 
star shows highly processed material at the surface (Paper I and Walborn et al. 2004), and 
in Paper I we argue that this star might be the result of binary evolution. We note that this 
star is found to the left of the ZAMS in Fig. 26, a nonphysical location. A similar problem 
was found for some of the relatively close binaries in the R136 cluster by Massey et al. (2002), 
suggesting that the stars had suffered from some binary interactions. However, we are also 
left with the possibility that at high Tes the models may be underestimating either logg or 
R] given that R is derived from My and Tes, a problem with logg appears to be the most 
likely. 


6. Summary, Discussion, and Future Work 

We have analyzed 40 0-type stars in the Magellanic Clouds, including many stars of 
the earliest type. Modeling was successful for 33 of these stars, with the other 7 showing 
the spectroscopic signature of unresolved companions. This study, in combination with an 
analysis of 24 Galactic stars by Repolust et al. (2004), using similar techniques and the same 
model atmosphere code, allow us to obtain the following results: 

1. The effective temperatures of 03-07 dwarfs and giants in the SMC is about 4000 K 
hotter than for stars of the same spectral type in the Milky Way. The differences 
decrease as one approaches BO V. This is readily understood in terms of the decreased 
importance of line- and wind-blanketing at the lower metallicity that characterizes the 
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SMC. The results for the LMC appear to be intermediate between the two galaxies. 
A similar effect is seen for the supergiants, although the differences decrease more 
rapidly with increasing spectral type (i.e., Table 8 and Fig. 27). For each galaxy, 
there is no difference in the effective temperatures of dwarfs and giants of the same 
spectral type, while the supergiants are about 4000 K cooler than the dwarfs for the 
hottest types in the SMC, and about 6000K cooler than the dwarfs for the hottest 
types in the Milky Way. The differences in effective temperatures between supergiants 
and dwarfs decrease for the later 0-types stars. Our effective temperature scale for 
dwarfs is in accord with some CMFGEN studies but not for others, and is signihcantly 
hotter than that indicated by WM-Basic modeling the UV. It is also hotter than the 
TLUSTY/CMFGEN modeling of NGC 346 stars by Bouret at al. (2003) and Heap et 
ah (2005). However, our scale leads to more consistent results with stellar evolutionary 
theory as evidenced by comparisons of the degree of coevality of stars in the NGC 346 
cluster (i.e. Fig. 30). 

2. Our data suggest that the wind momentum (T^mom = VoqMR^'^) scales with luminosity 
in the expected way with metallicity as predicted by radiatively driven wind theory 
(Kudritzki 2002, Vink et al. 2001): 

n rsr rl/“eff 

■^mom ^ 1 

with tteff ~ 0.55 fairly insensitive to Tgfj or Z, and with the constant of proportionality 
scaling with Two of the SMC stars with low mass-loss rates does not fit this 

relationship well, but most of the others do; a detailed analysis of the UV spectra (to 
be done in the future) might provide better constraints on the wind momenta in those 
cases. 

3. Most of the stars in our sample show a reasonable match between the spectroscopic 
mass and the evolutionary mass. However, stars with > 45000K show a systematic 
difference, with the spectroscopic mass signihcantly less (by a factor of 2 or more) 
than the evolutionary mass. This is similar to the long-standing “mass discrepancy” 
discussed by Herrero et al. (1992) for Galactic supergiants, but which has now been 
mostly resolved due to the lower effective temperatures of the Galactic models with 
the effects of line-blanketing included. Use of newer evolutionary tracks (which contain 
improved opacities, better treatment of mass-loss, and the effects of rotation) will tend 
to decrease the discrepancy, but such improvement is likely to be only of order 10%, 
and will not account for the factors of 2 discrepancies. We are left with the conclusion 
that the surface gravities or stellar radii may be underestimated in our models for stars 
of the highest effective temperatures. 
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4. We find that there is little correlation in the physical properties (such as Teg) with 
the new spectral types 02-3.5. Although stars in this group contain the hottest stars, 
neither the NIII/NIV emission line ratio nor even the He I/He H absorption line ratio 
provides a clear indication of the star’s effective temperature. For these extreme spec¬ 
tral types such line ratios appear to be sensitive to loggf and M as well as Teg. Thus 
although we can predict a star’s effective temperature from its two-dimension spectral 
type and metallicity for 04’s and later, it requires a detailed analysis of the entire 
spectrum to derive a reliable effective temperature (i.e., the mass-loss rate inferred 
from Ha and the terminal velocity of the wind from UV measurements are needed in 
addition to the blue optical spectrum. 

5. Two of the stars with the largest radial velocities with respect to the LMC are 03 stars 
listed by Massey et al. (1995) as “held” stars i.e., they are found far from the nearest 
OB association. These provide the hrst compelling evidence for the origin of this held 
population. 

Here we wish to emphasize again the point raised in the introduction that although we 
have used the best available data and (we believe) the best models and methodology, that 
the absolute ehective temperature scale of O stars is likely to undergo revision in the future 
as the physics is improved. Indeed, our study here underscores that other “best” methods, 
such as htting the UV lines via WM-BASIC, yield temperatures for Galactic stars that 
are signihcant cooler than that found by FASTWIND and some CMFGEN modeling of the 
optical spectra of an admittedly different sample. Reality is that different lines (and different 
spectral domains) seem to give different answers at this time, very likely independent of the 
model code used. That tells us that the physics used in our models is still imperfect. It 
is hard to tell which method is closer to the truth. Both have their merits. However, the 
strength of the current study is that we have been able to do a strictly differential comparison 
based on one and the same method applied to a large sample of objects in three galaxies of 
different metallicity. While the absolute scale might still be uncertain, we believe that the 
relative differences of the scales at the different metallicities obtained in our approach are 
accurate and signihcant. 

Future work should attempt to identify where improvements are needed; the logical 
place to start is by modeling the UV spectra of the stars studied here, and by modeling 
optical data on samples that have been analyzed only in the UV. In some cases we believe 
that better data of the same stars (such as the NGG 346 sample studied by Bouret et al. 
2003) may help resolve differences seen between our effective temperature scale and others, 
in much the same way as our better (nebular-subtracted) data on the R136 stars led to much 
better model hts (Paper I). We keenly anticipate the results ahead. 
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Table 1. Program Stars®" 


Name^ 

Cat ID‘= 

0:2000 

52000 

V 

to 

1 

U-B 

E{B - 

Mv-® 

Spectral Type^ 

AV 177 

SMC-038024 

00 56 44.17 

-72 03 31.3 

14.53 

-0.21 

-1.05 

0.12 

-4.78 

04 V((f)) 

AV 435 

SMC-067670 

01 08 17.88 

-71 59 54.3 

14.00 

-0.06 

-0.98 

0.28 

-5.81 

03 V((f*)) 

AV 440 

SMC-068756 

01 08 56.01 

-71 52 46.5 

14.48 

-0.18 

-1.00 

0.15 

-4.93 

08 V 

AV 446 

SMC-069555 

01 09 25.46 

-73 09 29.7 

14.59 

-0.24 

-1.06 

0.10 

-4.66 

06.5 V 

AV 476 

SMC-074608 

01 13 42.41 

-73 17 29.3 

13.52 

-0.09 

-0.93 

0.28 

-6.29 

02-3 V -1- comp. 

Sk-67°22=BAT99-12 

LMC-034056 

04 57 27.47 

-67 39 03.3 

13.44 

-0.18 

-1.05 

0.16 

-5.53 

02 If* 

Sk-65°47=LH43-18 


05 20 54.67 

-65 27 18.3 

12.68 

-0.13 

-0.93: 

0.19 

-6.39 

04 If 

LH58-496=LH58-10aS 


05 26 44.21 

-68 48 42.1 

13.73 

-0.23 

-1.09 

0.11 

-5.09 

05 V(f) 

LH81:W28-23 


05 34 50.11 

-69 46 32.3 

13.81 

-0.16 

-1.13 

0.15 

-5.14 

03.5 V((f-k)) 

LH90: Br58=B AT99-68‘' 


05 35 42.02 

-69 11 54.2 

14.13 

0.49 

-0.48 

0.85 

-6.98 

03If/WN6 

LH90:ST2-22 


05 35 45.26 

-69 11 35.1 

13.93‘ 

0.18' 

-0.73' 

0.60 

-6.41 

03.5 Ill(f-I-) 

BI 237 

LMC-164942 

05 36 14.68 

-67 39 19.3 

13.89 

-0.12 

-0.97 

0.25 

-5.38 

02 V((f*)) 

BI 253 

LMC-168644 

05 37 34.48 

-69 01 10.4 

13.76 

-0.09 

-1.02 

0.25 

-5.50 

02 V((f*)) 

LH101:W3-14=ST5-52-i 


05 39 05.41 

-69 29 20.7 

13.41 

-0.15 

-0.89 

0.17 

-5.60 

03 V -1- 0 V Composite 

LH101:W3-19=ST5-31 


05 39 12.20 

-69 30 37.6 

12.37‘ 

-0.06' 

-0.92' 

0.30 

-7.04 

02 If* 

LMC2-675 


05 43 13.00 

-67 51 16.0 

13.66 

-0.25 

-1.14 

0.07 

-5.04 

02 III(f*) -1- 0 V Composite 

R136-007=Mk39'^ 


05 38 40.3186 

-69 06 00.172 

13.01 



0.47 

-6.9 

02 If* Composite 

R136-014=Mk37'^ 


05 38 42.4986 

-69 06 15.396 

13.57 



0.48 

-6.4 

03.5 If* 

R136-018'' 


05 38 44.2211 

-69 05 56.954 

13.91 



0.42 

-5.9 

03 III 

R136-033'' 


05 38 42.2106 

-69 06 00.988 

14.43 



0.35 

-5.1 

03 V 


^Coordinates and photometry are from Massey 2002 or Massey, Waterhouse, Sz DeGioia-Eastwood 2002 unless otherwise noted. 

'^Identifications are as follows: “AV” from Azzopardi Sz, Vigneau 1982; “BI” from Brunet et al. 1975; “BAT” from Breysacher, Azzopardi, Sz Testor 1999; 
“Br” from Breysacher 1981; “LH” from Lucke 1972 except as noted, “LMC2” from Massey et al. 1995 and Massey 2002; “Sk” from Sanduleak 1969; “ST” 

from Testor & Niemela 1998; “W” from Westerlund 1961; “Mk” from Melnick 1985; “R136-NNN” from Hunter et al. 1997 and Massey & Hunter 1998 

^Designations from the catalog of Massey 2002. 

^From averaging the color excesses in B — V and U — B based upon the spectral type. See Massey 1998b. 

^Computed using Ay = 3.1 x E(B — V), with assumed distance moduli for the SMC and LMC of 18.9 and 18.5, respectively (Westerlund 1997, van den 

Bergh 2000). 

^New to this paper. 

SH58-10A identification is from Lucke 1972; LH58-496 identification is from Garmany at al. 1994. 

'^Misidentified with a neighboring bright star in Table 4 of Massey 2002. 













^Photometry new to this paper, based upon the CCD material described by Massey 2002. 
•^Coordinates and photometry from Testor Sz Niemela 1998. 

^Coordinates and photometry from Hunter et al. 1997 and Massey Sz Hunter 1998. 
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Table 2. Sources of Data Used in This Study 


Spectral Region 

Telescope 

Instrument 

Aperture 
(arcsec x arcsec) 

Grating 

Wavelength 

(A) 

Resolution 

(A) 

S/N^ 

Use 

Far-UV 

FUSE/C002 


30x30 


905-1187 

0.1 

40 

Terminal velocities 

UV 

HST/9412 

STIS 

0.2x0.2 

G140L 

1150-1740 

0.9 

30 

Terminal velocities 

Near-UV 

HST/9412 

STIS 

0.2x52 

G430M 

3020-3300 

0.4 

10-80 

Check on the modeling 

Blue-optical 

CTIO 4-m 

RC Spec 

1.3x330 

KPGLD 

3750-4900 

1.4 

400-500 

Modeling 


HST/7739,9412 

STIS/CCD 

0.2x52 

G430M 

4310-4590 

0.4 

100 

Modeling 


HST/6417 

FOS 

0.26(circ) 

G400M 

3250-4820 

3.0 

40-60 

Spectral class R136 stars 

Ho: 

CTIO 4-m 

RC Spec 

1.3x330 

KPGLD 

5400-7800 

2.8 

150 

Modeling 


HST/7739,9412 

STIS/CCD 

0.2x52 

G750M 

6300-6850 

0.8 

60 

Modeling 


^Signal to noise per spectral resolution element. 




Table 3. 

Spectral Regions Observed 


Star 

FUV 

UV 

NUV 

Blue-optical 

Ho 

AV 177 

FUSE 

HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

AV 435 


HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

AV 440 


HST/9412 


CTIO 4-m 

HST/9412 

AV 446 


HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

AV 476 

FUSE 

HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

Sk-67°22 

FUSE 

HST/9412 


CTIO 4-m 

HST/9412 

Sk-65°47 

FUSE 

HST/9412 

HST/9412 

CTIO 4-m, HST/9412 

HST/9412 

LH58-496 


HST/9412 

HST/9412 

CTIO 4-m 

HST/9412 

LH81:W28-23 


HST/9412 

HST/9412 

CTIO 4-m 

CTIO 4-m, HST/9412 

LH90:Br58 


HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

LH90:ST2-22 


HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

BI 237 


HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

BI 253 

FUSE 

HST/9412 


CTIO 4-m 

CTIO 4-m, HST/9412 

LH101:W3-14 


HST/9412 

HST/9412 

CTIO 4-m 

CTIO 4-m, HST/9412 

LH101:W3-19 


HST/9412 

HST/9412 

CTIO 4-m, HST/9412 

CTIO 4-m, HST/9412 

LMC2-675 


HST/9412 

HST/9412 

CTIO 4-m 

CTIO 4-m, HST/9412 

R136-007 


HST/9412 


HST/7739, HST/6417 

HST/7739 

R136-014 


HST/9412 


HST/7739, HST/6417 

HST/7739 

R136-018 


HST/9412 


HST/7739, HST/6417 

HST/7739 

R136-033 


HST/9412 


HST/7739, HST/6417 

HST/7739 













Table 4. Terminal Velocities in km s ^ 


Star 

'^OO 

Lines Used 

Comments 

AV 177 

2650 

CIV, NV 


AV 435 

1500:: 

CIV, SilV 

Weak CIV 

AV 440 

1300: 

CIV 

Weak CIV 

AV 446 

1400: 

CIV 

Weak CIV 

AV 476 

2670 

CIV 


Sk-67°22 

2650 

CIV 


Sk-65°47 

2100 

CIV 


LH58-496 

2400:: 

CIV 

Wide CIV, plus abs. contamination 

LH81:W28-23 

3050 

CIV 


LH90:Br58 

1900: 

SilV 

Could not use CIV 

LH90:ST2-22 

2560 

CIV 


BI 237 

3400 

CIV 


BI 253 

3180 

CIV 


LH101:W3-14 

3100 

CIV 


LH101:W3-19 

2850 

CIV 


LMC2-675 

3200 

CIV 


R136-007 

2100 

CIV, SilV 


R136-014 

2000 

CIV, SilV 

Velocity scale uncertain 

R136-018 

3200 

CIV 


R136-033 

3250 

CIV 







Table 5. Results of Model Fits 


Name 

Spectral 

Type 

Ustf 

(°K) 

log 9eff 
[egs] 

Strue^ 

[cgs] 

R 

(Rq) 

Mv 

mags 

BC 

mags 

Mbol 

mags 

Mass 

M© 

M 

(IO-^Mq yr-l) 

/3 

Uoo 

(km s-1) 

He/H'= 

Comments 

AV 177 

04 V((f)) 

44000 

3.80 

3.85 

8.9 

-4.78 

-4.04 

-8.82 

21 

0.3 

0.8 

2650 

0.15 

M < 0.5 

AV 435 

03 V((f*)) 

45000 

3.80 

3.81 

14.2 

-5.81 

-4.12 

-9.93 

48 

0.5 

0.8 

1500:: 

0.10 


AV 440 

08 V 

37000 

4.00 

4.01 

10.6 

-4.93 

-3.52 

-8.45 

42 

0.1 

0.8 

1300: 

0.12 

M < 0.3 

AV 446 

06.5 V 

41000 

4.15 

4.15 

00 

bo 

-4.66 

-3.82 

-8.48 

40 

0.1 

0.8 

1400 

0.15 

M < 0.3 

AV 476 

02-3 V-|-comp. 




-6.29 







2670 


Composite 

Sk -67°22 

02 If* 

> 42000 

«3.5 

«3.56 

SS13.2 

-5.53 

< -3.94 

-9.47 

^23 

15 

0.8 

2650 

0.30 

Tgff lower limit 

Sk -65°47 

04 If 

40000 

3.60 

3.62 

20.1 

-6.39 

-3.79 

-10.18 

61 

12 

0.8 

2100 

0.10 


LH58-496 

05 V 

42000 

4.00 

4.04 

10.5 

-5.09 

-3.88 

-8.97 

44 

0.6 

0.8 

2400 

0.10 


LH81:W28-23 

03.5 V((f+)) 

47500 

3.80 

3.81 

10.0 

-5.14 

-4.26 

-9.40 

24 

2.5 

0.8 

3050 

0.25 


LH90:Br58 

03If/WN6 

40000-42000: 

3.5: 

3.5: 

30: 

-6.98 

-4.1:: 

-11.1: 

100:: 

40: 

3.0 

1900 

0.1:: 

See text 

LH90:ST2-22 

03.5 III(f+) 

44000 

3.70 

3.71 

18.9 

-6.41 

-4.05 

-10.46 

67 

4.5 

0.8 

2560 

0.20 


BI 237 

02 V((f*)) 

48000 

3.90 

3.92 

11.1 

-5.38 

-4.29 

-9.67 

37 

2.0 

0.8 

3400 

0.10 


BI 253 

02 V((f*)) 

> 48000 

^3.90 

«3.93 

^11.8 

-5.50 

<-4.30 

<-9.80 

^43 

3.5 

0.8 

3180 

0.10 

Tgff lower limit 

LH101:W3-14 

03 V + O V 




-5.60 







3100 


Composite 

LH101:W3-19 

02 If* 

> 44000 

^3.90 

«3.91 

^ 25.5 

-7.04 

<-4.06 

<-11.10 

^193 

20 

0.8 

2850 

0.10 

Tgff lower limit 

LMC2-675 

02 III(f*) + O V 




-5.04 







3200 


Composite 

R136-007 

02 If* Composite 




-6.9 







2100 


Composite 

R136-014 

03.5 If* 

38000 

3.50 

3.51 

21.1 

-6.4 

-3.65 

-10.0 

53 

23 

0.8 

2000 

0.10 


R136-018 

03 III{f*) 

45000 

3.75 

3.77 

14.7 

-5.9 

-4.11 

-10.0 

46 

2.0 

0.8 

3200 

0.10 


R136-033 

02-3.5 V 

47000 

3.75 

3.77 

9.8 

-5.1 

-4.23 

-9.3 

21 

2.0 

0.8 

3250 

0.10 

1 




By number. 
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Table 6. OVI vs. CIV Terminal Velocities (km s 


Star 

Sp.Type 

OVI 

CIV 

AV 177 

04 V((f)) 

2500 

2650 

Sk-67°22 

02 If* 

2650 

2650 

Sk-65°47 

04 If 

2000 

2100 

BI 253 

02 V((f*)) 

3200 

3180 






Table 7. Line Strengths in the Earliest 0-type Stars 


Star 

Type 

^eff 

(1000 K) 

S'true 

(cgs) 

M 

(lO-^M© yr-l) 

EWs 

He I A4471 

[mA]“ 

He II A4542 

, EW(Hel) 

log EW(Rell) 

EWs 

N IV A4058 

[mA]'= 

N III A4634, 42 

, EW(NIII) 

log E:W(N1V) 

Comments 

DWARFS 

BI 237 

02 V((f*)) 

48.0 

3.92 

2.0 

20 ± 10 

750 ± 

10 

-1.6 ± 0.2 

-165 ± 20 

-90 ± 20 

-0.3 ± 0.1 

C IV A 4658 

BI 253 

02 V((f*)) 

>48.0 

3.93 

3.5 

< 25 ± 10 

710 ± 

10 

< -1.5 ± 0.2 

-450 ± 20 

> -60 ± 20 

< -0.9 ± 0.1 

C IV A 4658 

AV 435“ 

03 V((f*)) 

45.0 

3.81 

0.5 

125 ± 10 

750 ± 

10 

-0.8 ± 0.0 

-80 ± 20 

> -60 ± 20 

< -0.1 ± 0.2 


LH101:W3-24 

03.5 V((f+)) 

48.0 

4.01 

0.5 

120 ± 10 

790 ± 

10 

-0.8 ± 0.0 

> -20 ± 20 

-200 ± 20 

> 1.0 ± 0.4 

03 V((f)) in Pap I. 

LH81:W28-23 

03.5 V((f+)) 

47.5 

3.81 

2.5 

80 ± 10 

830 ± 

10 

-1.0 ± 0.1 

-120 ± 20 

-500 ± 20 

0.6 ± 0.1 


AV 177“ 

04 V((f)) 

44.0 

3.85 

0.3 

170 ± 10 

950 ± 

10 

-0.8 ± 0.0 

> -50 ± 20 

-250 ± 20 

> 0.7 ± 0.2 


LH81:W28-5 

04 V((f+)) 

46.0 

3.81 

1.2 

200 ± 10 

890 ± 

10 

-0.7 ± 0.0 

-50 ± 20 

-480 ± 20 

1.0 ± 0.2 


R136-055 

02-3.5 V 

47.5 

3.81 

0.9 

35 ± 20 

650 ± 

20 

-1.3 ± 0.2 

> -200 ± 100 

> -200 ± 100 


03 V in Paper I 

R136-033 

02-3.5 V 

47.0 

3.77 

2.0 

50 ± 20 

790 ± 

20 

< -1.2 ± 0.2 

> -200 ± 100 

> -200 ± 100 


03 V in this paper 

R136-040 

02-3.5 V 

>51.0 

^3.81 

2.0 

< 100 ± 50 

760 ± 

20 

< -0.9 ± 0.2 

> -200 ± 100 

> -200 ± 100 


03 V in Paper I 

GIANTS 

LH64-16 

ON2 III(P) 

54.5 

3.91 

4.0 

100 ± 10 

1090 ± 10 

-1.0 ± 0.0 

-660 ± 20 

-120 ± 20 

-0.7 ± 0.1 

C III A 46^8 

R136-047 

02 III(f*) 

>51.0 

^3.91 

6.0 

< 30 ± 20 

450 ± 

20 

< -1.2 ± 0.3 

-700 ± 100 

> -250 ± 100 

< -0.5 ± 0.2 


R136-018 

03 Ill(f^) 

45.0 

3.77 

2.0 

90 ± 20 

670 ± 

20 

-0.9 ± 0.1 

-750 ± 100 

> -150 ± 100 

< -0.7 ± 0.3 

CD 

LH90:ST2-22 

03.5 III(f+) 

44.0 

3.71 

4.5 

140 ± 10 

00 

00 

o 

10 

-0.9 ± 0.0 

-230 ± 20 

-750 ± 20 

0.5 ± 0.0 


SUPERGIANTS 

Sk-67 22 

02 If* 

>42.0 

SS3.56 

15 

< 50 ± 20 

650 ± 

20 

< -1.1 ± 0.2 

-1610 ± 20 

-350 ± 20 

-0.7 ± 0.0 

C IV A4658 

LH101:W3-19 

02 If* 

>44.0 

^3.91 

<20 

<- 50 ± 20 

570 ± 

20 

< -1.1 ± 0.2 

-300 ± 20 

-230 ± 20 

-0.1 ± 0.1 

C IV A4658 

R136-020 

02 If* 

>42.5 

^3.61 

23 

< 30 ± 20 

510 ± 

20 

< -1.2 ± 0.3 

-2570 ± 50 

> -100 ± 50 

< -1.4 ± 0.2 


R136-036 

02 If* 

>43.0 

^3.71 

14 

< 50 ± 20 

550 ± 

20 

< -1.0 ± 0.2 

-2020 ± 100 

> -200 ± 100 

< -1.0 ± 0.2 


LH90:Br58'’ 

03 If/WN6 

40-42 

3.5: 

40: 

< 50 ± 20 

o 

o 

20 

< -0.9 ± 0.2 

-1010 ± 20 

-480 ± 20 

-0.3 ± 0.0 


R136-014 

03.5 If* 

38.0 

3.51 

23 

40 ± 20 

o 

o 

20 

-1.0 ± 0.2 

-1360 ± 100 

-1400 ± 100 

0.0 ± 0.0 


Sk-65° 47 

04 If 

40.0 

3.62 

12 

115 ± 10 

550 ± 

10 

-0.7 ± 0.0 

-175 ± 20 

-1540 ± 20 

0.9 ± 0.1 



^SMC star. 

^This star is excluded from consideration in Fig. 24 due to the uncertainity in its effective temperature. 







Table 8. Effective Temperature Scale [K] 


Type 

Milky Way 

SMC 

V+III 

I 

V+III 

I 

03 

46500 

40250 

49500 

45250 

04 

44000 

39000 

47750 

43000 

05 

41000 

37750 

45000 

41000 

05.5 

39500 

36750 

43500 

40000 

06 

38250 

36000 

42250 

38500 

06.5 

37000 

35500 

41000 

37500 

07 

36000 

34750 

39250 

36250 

07.5 

34750 

34000 

37750 

35250 

08 

33750 

33000 

36250 

34000 

08.5 

32750 

32500 

34500 

33000 

09 

31750 

31750 

33000 

32000 

09.5 

30750 

30750 

31500 

30750 

BO 

30000 

29750 

30000 

29750 







Table 9. Parameters for the Complete Sample 


Star 

Spectral 

Type 

V sin i 
km s~ ^ 

T"eff 

[1000 K] 

5true 

[cgs] 

M[10“'= 

yr“^] 

Uoo 

[km s“^] 

R/Rq 

He/H 

[num.] 

log Dj-nom 

log L/Lq 

•Mgpect 

-^evol 

DWARFS 

R136-040 

02-3.5 V 

120 

>51.0 

«3.81 

2.0 

3400 

sslO.3 

0.10 

?s29.14 

>5.82 

ss25 

>70 

BI 253 

02 V((f*)) 

200 

>48.0 

^3.93 

3.5 

3180 

^11.8 

0.10 

?s29.38 

>5.82 

ss43 

>64 

BI 237 

02 V((f*)) 

150 

48.0 

3.92 

2.0 

3400 

11.1 

0.10 

29.16 

5.77 

37 

62 

R136-033 

02-3.5 V 

120 

47.0 

3.77 

2.0 

3250 

9.8 

0.10 

29.11 

5.6 

21 

52 

R136-055 

02-3.5 V 

120 

47.5 

3.81 

0.9 

3250 

9.4 

0.10 

28.75 

5.6 

21 

52 

LH101:W3-24 

03 V((f)) 

120 

48.0 

4.01 

0.5 

2400 

8.1 

0.15 

28.33 

5.49 

25 

48 

AV 435'^ 

03 V((f*)) 

110 

45.0 

3.81 

0.5 

1500:: 

14.2 

0.10 

28.3:: 

5.87 

48 

62 

LH81:W28-23 

03 V((f+)) 

120 

47.5 

3.81 

2.5 

3050 

10.0 

0.20 

29.18 

5.66 

24 

55 

AV 177'= 

04 V((f)) 

220 

44.0 

3.85 

0.3 

2650 

8.9 

0.15 

28.18 

5.43 

21 

40 

LH81:W28-5 

04 V((f+)) 

120 

46.0 

3.81 

1.2 

2700 

9.6 

0.20 

28.80 

5.58 

22 

50 

AV 377'= 

05 V((f)) 

120 

45.5 

4.01 

<0.3 

2350 

9.1 

0.35 

<28.13 

5.52 

31 

45 

AV 14'= 

05 V 

150 

44.0 

4.01 

<0.3 

2000 

14.2 

0.10 

<28.15 

5.85 

75 

59 

LH58-496 

05 V 

250 

42.0 

4.04 

0.6 

2400 

10.5 

0.10 

28.47 

5.49 

44 

41 

AV 446'= 

06.5 V 

95 

41.0 

4.15 

<0.3 

1400 

00 

bo 

0.15 

<27.90 

5.29 

40 

33 

AV 207'= 

07.5 V((f)) 

120 

37.0 

3.72 

<0.3 

2000 

11.0 

0.10 

<28.10 

5.32 

23 

29 

AV 296'= 

07.5 V((f)) 

300 

35.0 

3.63 

0.5 

2000 

11.9 

0.10 

28.34 

5.28 

22 

29 

AV 440'= 

08 V 

100 

37.0 

4.01 

<0.3 

1300: 

10.6 

0.12 

<27.9: 

5.28 

42 

30 

GIANTS 

LH64-16 

ON2 III(P) 

120 

54.5 

3.91 

4.0 

3250 

9.4 

1.0 

29.40 

5.85 

26 

76 

R136-047 

02 III(f*) 

120 

>51.0 

^3.91 

6.0 

3500 

sslO.4 

0.10 

^28.64 

>5.82 

:^32 

>70 

R136-018 

03 III(P) 

180 

45.0 

3.77 

2.0 

3200 

14.7 

0.10 

29.19 

5.9 

46 

65 

LH90:ST2-22 

03.5 III(f+) 

120 

44.0 

3.71 

4.5 

2560 

18.9 

0.20 

29.50 

6.08 

67 

79 

AV 378'= 

09.5 III 

110 

31.5 

3.27 



15.4 

0.15 


5.34 

17 

29 

AV 396'= 

BO III 

120 

30.0 

3.52 



14.1 

0.15 


5.17 

24 

22 

SUPERGIANTS 

LH101:W3-19 

02 If* 

180 

>44.0 

«3.91 

20' 

2850 

^25.5 

0.10 

JS30.26 

>6.34 

=2193 

>109 

R136-036 

02 If* 

120 

>43.0 

==;3.71 

14° 

3700 

^12.8 

0.10 

?=;30.07 

>5.7 

:^31 

>50 

R136-020 

02 If* 

120 

>42.5 

^3.61 

23° 

3400 

^16A 

0.20 

?=;30.30 

>5.9 

>40 

>62 

Sk -67 22 

02 If* 

200 

>42.0 

^3.56 

15^ 

2650 

^13.2 

0.30 

^29.96 

^5.69 

>23 

>49 

LH90:Br58 

03If/WN6 


40.0-42.0 

3.5: 

40:° 

1900 

30: 

30.4 

0.1:: 

6.3: 

40: 

101: 

R136-014 

03.5 If* 

120 

38.0 

3.51 

23° 

2000 

21.1 

0.10 

30.12 

5.9 

53 

57 

Sk -65 47 

04 If 

160 

40.0 

3.62 

12° 

2100 

20.1 

0.10 

29.85 

5.97 

61 

65 

AV 75'= 

05.5 1(f) 

120 

40.0 

3.61 

3.5 

2100 

25.4 

0.10 

29.37 

6.16 

96 

84 

AV 26'= 

06 1(f) 

150 

38.0 

3.52 

2.5 

2150 

27.5 

0.10 

29.25 

6.14 

91 

81 

AV 469'= 

08.5 1(f) 

120 

32.0 

3.13 

1.8 

2000 

21.2 

0.20 

29.02 

5.64 

22 

39 


''From the non-rotating models of Charbonnel et al. 1993 (SMC) and Schaerer et al. 1993 (LMC). 

‘SMC 

Ho: in emission; these values for M should be reduced by a factor of 0.44 to allow for the effects of clumping in the stellar wind. See text. 
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'n'QveleiL^h (Angstroms) 



Fig. 1.— AV 177. (a) A portion of the blue-optical spectrum of AV 177 is shown with the 
major lines identihed. (b) Selected spectral lines (black) are shown compared to the model 
fits (red in the on-line version; thin black line in the printed version). The bar to the left of 
each line shows a change of 20% intensity relative to the continuum, and the top of the bar 
denotes the continuum level. A radial velocity of 100 km s“^ and a rotational broadening 
nsini of 220 km s“^ was used in making this comparison. 
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'CTavelen^b (Angstroms) 



Fig. 2.— AV 435. Same as Fig. 1. A radial velocity of 200 km s ^ and a rotational broadening 
nsini of 110 km was nsed in making this comparison. 
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Fig. 3.— AV 440. Same as Fig. 1. A radial velocity of 190 km s ^ and a rotational broadening 
nsini of 100 km was nsed in making this comparison. 
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Fig. 4.— AV 446. Same as Fig. 1. A radial velocity of 140 km s ^ and a rotational broadening 
nsini of 95 km was nsed in making this comparison. 
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Fig. 6.— Sk—67°22. Same as Fig. 1. A radial velocity of 430 km and a rotational broad¬ 
ening nsini of 200 km was nsed in making this comparison. The HST Ha observation 
is shown. 
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Fig. 7.— Sk—65°47. Same as Fig. 1. A radial velocity of 250 km s ^ and a rotational 
broadening nsini of 160 km was nsed in making this comparison. 
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'n’avelen^b (Angstroms) 



Fig. 8.— LH58-496. Same as Fig. 1. A radial velocity of 275 km s ^ and a rotational 
broadening nsini of 250 km was nsed in making this comparison. 
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'n'Qveleii^b (Angstroms) 



Fig. 9.— LH81:W28-23. Same as Fig. 1. A radial velocity of 350 km and a rotational 
broadening nsini of 120 km was used in making this comparison. A cosmic ray on the 
wing of the He II A4200 has been removed in this plot. 
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Fig. 10.— LH90:Br58. The strong emission in this star prevented a satisfactory ht, although 
tentative parameters are given in the text. Note the strong diffuse interstellar band at A4430. 
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'n'Qveleii^b (Angstroms) 



Fig. 11.— LH90:ST2-22. Same as Fig. 1. A radial velocity of 250 km and a rotational 
broadening nsini of 120 km was used in making this comparison. A cosmic ray on the 
He II A4200 prohle has been removed in this £t. 
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'n'Qveleii^b (Angstroms) 



Fig. 12.— BI 237. Same as Fig. 1. A radial velocity of 430 km and a rotational broadening 
nsini of 150 km was used in making this comparison. The He II A4200 prohle was is 
flat-bottomed, presumably due to a reduction problem, and not used in making the £t. 
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Fig. 13.— BI 253. Same as Fig. 1. A radial velocity of 270 km s ^ and 
broadening nsini of 200 km was nsed in making this comparison. 
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Fig. 14.— LH101:W3-14. The spectrum shown here is composite, with an 03 star dominat¬ 
ing the nitrogen and He II spectra, and a later-type O star contributing the He I. No single 
model was able to £t the strengths of the He I and He H lines. 
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'n’avelen^b (Angstroms) 



Fig. 15.— LH101:W3-19. Same as Fig. 1. A radial velocity of 320 km s ^ and a rotational 
broadening nsini of 180 km was used in making this comparison. 
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Fig. 16.— LMC2-675. The spectrum shown here is composite, with an 03 star dominating 
the nitrogen and He II spectra, and a later 0-type star contributing the He I. As in Fig. 14, 
no single model was able to £t the strengths of the He I and He H lines. 





Fig. 17.— R136-007. We identify the major lines in the spectrum of R136-007. The data 
are from the FOS observations of Massey & Hunter (1998), except for the region 4310A 
to 4590A, where we have spliced in our higher S/N STIS spectrum. No combination of 
parameters led to a good £t, and we conclude that the star is composite, consistent with the 
discovery of eclipses in the light-curve by Massey et al. (2002). 
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Fig. 18.— R136-014. The data in (a) come from the FOS observations of Massey & Hunter 
(1998), except for the region 4310A to 4590A, where we have spliced in our higher S/N STIS 
spectrum. In (b) the prohles are from the STIS data, except for that of He H A4686, which 
comes from the FOS data. The colors and symbols have the same meaning as in Fig. 1. Prob¬ 
lems with the STIS wavelength zero-point preclude determining an accurate radial velocity; 
a rotational broadening nsini of 120 km s“^ was used in making this comparison. 
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Fig. 19.— R136-018. Same as Fig. 18. A rotational broadening nsini of 180 km s ^ was 
used in making this comparison. 
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'n’avelen^b (Angstroms) 



Fig. 20.— R136-033. Same as Fig. 18. A rotational broadening nsini of 120 km s ^ was 
used in making this comparison. 
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Fig. 21.— CIV vs OVI Terminal Velocities. For the star Sk—67°22 we show CIV prohle 
(upper panel) and OVI prohle (lower panel) along with our adopted SEI hts. In both cases 
we found a terminal velocity v^o = 2650 km s“^. The dotted curves (blue in the on-line 
version) show the sensitivity of the velocity to a change of ±100 km s“^. Since there has 
been no correction for underlying photospheric absorption for CIV, the poor £t on the long 
wavelength side of CIV is expected; it is the short wavelength side which dehnes the terminal 
velocity. Similarly the agreement with the OVI could be improved by better correction for 
interstellar absorption, but the short wavelength side demonstrates that the terminal velocity 
is well determined. Note the presence of Lyman (3 emission in the OVI profile. 
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Wavelength (A) 


Fig. 22.— The NUV region. The region containing the He I A3187 and He H A3203 lines are 
shown for the four non-composite stars in our sample for which we have data in this region. 
Although we did not use this region in obtaining the hnal parameters, each star shows good 
agreement with the corresponding model. 
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Fig. 23.— NIII/NIV vs Hel/Hell for 02-3.5 stars. The log of the ratio of the equivalent 
widths (EWs) of N III A4634,42 to N IV A4058 emission is shown as a function of the log of 
the ratio of the EWs of He I A4471 to He H A4542. Typical errors (as given in Table 7) are 
0.2 dex on each axis, but are not shown to prevent confusion with lower and upper limits. 
Dwarfs are shown as hlled circles, giants are shown as hlled triangles, and supergiants are 
shown as open circles. Red symbols are for the LMC stars; green for the SMC. 
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Fig. 24.— Effective temperatures of the earliest-type O stars (02-03.5) is shown as a 
function of (a) the Hel/Hell line ratio, and (b) the NIII/NIV emission line ratio. Again, the 
uncertainty in the line ratios are typically 0.2 dex; see Table 7. The symbols have the same 
meaning as in Fig. 23. 
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Fig. 25.— The dependence of loghF' = log[EW(He IA4471)/(He IIA4542) on effective tem¬ 
perature and mass-loss rates according to FASTWIND. The contours show values of constant 
loghF' (green for positive, black for negative, with the values indicated). In (a) we show 
the dependence of log W' on log g and Tgs with the value of log Q (as defined below) held 
constant (-12.45, correspond to a mass-loss rate of yr“^ for Vca = 2000 km s“^ for a 

star with radius lOi?©). Although log W’ is primarily a temperature indicator over most of 
parameter space, it becomes a sensitive function of surface gravity for the hottest stars, (b) 
There is also a strong dependence of log W on the mass-loss rates for hot stars with high 
rates. The value Q is defined as M/{voo * RY'^■ The data here all have log^f = 3.8. 
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Fig. 26.— The location of our stars in the H-R diagram. Dwarfs are indicated by hlled 
circles, giants by triangles, and supergiants by open circles, (a) The non-rotating stellar 
evolutionary models of Charbonnel et ah (1993) are shown for the SMC, and (b) those of 
Schaerer et al. (1993) are shown for the LMC. Only the H-burning part of the tracks are 
shown. For comparison, we include the newer rotation models (dotted lines) for initial masses 
of 6 OM 0 and 4 OM 0 . 
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Fig. 27.— The effective temperatures as a function of spectral subtype, with a value of 2 
denoting spectral type 02, 5.5 denoting 05.5, and “10” denoting BO. The filled circles for 
dwarfs; the filled triangles are giants, and the open circles are supergiants. Black symbols 
correspond to the Milky Way, red to the LMC, and green to the SMC. The data for the 
Milky Way is taken from Repolust et al. (2004); the data for the SMC and LMC are taken 
from Paper I and the present study. 
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Fig. 28.— The effective temperature scale adopted here for the SMC (green) and the Milky 
Way (black). Solid lines denote the scale for supergiants; dashed lines are for dwarfs and 
giants, (a) The supergiant and dwarf effective temperature scales from Vacca et al. (1996) are 
shown for comparison in blue, (b) The effective temperature scale is shown with comparison 
for the data for the dwarfs and giants, (c) The supergiant effective temperature scale is 
shown in comparison with the data. The LMC scale is assumed to be intermediate between 
that of the SMC and the Milky Way. 
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Fig. 29.— Our effective temperature scales (dashed or solid curves) are compared to recent 
data from others. Black symbols represent the Milky Way; red, the LMC; and green, the 
SMC. The different symbols are indicated in the key, where the data come from Bianchi 
& Garcia (2002), [BG02]; Garcia & Bianchi (2004), [GB04]; Growther et ah (2002), [G02]; 
Hillier et ah (2003), [H03]; and Herrero et ah (2002), [HPN02], Bouret et al. (2003), [B03], 
Martins et al. (2002), [MSH02], Martins et al. (2004), [MSH04], and Evans et al. (2004) 
[E04]. 
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Fig. 30.— . The H-R diagram for the NGC 346 stars studied by Bouret et ah (2003), 
excluding NGC 346-MPG12, which Walborn et ah (2000) argue is not part of the cluster. 
The evolutionary tracks (solid lines) are from Charbonnel et al. (1993) computed for a 
metallicity of Z=0.004. The dotted lines are isochrones computed form these models for 
ages of 1-10 Myr at intervals of 1 Myr. (a) Stars placed with the parameters of Bouret et al. 
(2003). (b) Stars placed using the average spectral type to effective temperature calibration 
given here. In both cases the same photometry and color excesses were used, although a 
more realistic estimate of the reddening would result in stars about 0.1 mag more luminous. 
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Effective Temperature (1000 K) 


Fig. 31.— The bolometric corrections are shown as a fnnction of effective temperatnre. The 
black points are for the Milky Way; the red points are for the LMC, and the green points 
are for the SMC. The smooth cnrve denotes the relationship given in Eqn. 2. 
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Fig. 32.— The wind momentum luminosity relationship. The symbols have the same mean¬ 
ing as in Fig. 27. Stars with only lower limits on their effective temperatures are not shown, 
and stars whose values are particularly uncertain are indicated by parenthesis. Typical error 
bars are shown in the upper left, (a) The data are shown uncorrected for wind clumping, 
(b) The stars showing Ha emission have been corrected by -0.36 dex in T^mom to correct for 
the effects of wind clumping on the deduced mass-loss rates. The three lines are not fits to 
the data; instead, these are the theoretical expectations from radiatively driven wind theory. 
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Fig. 33.— The evolutionary masses are compared to the spectroscopic masses. The circles 
represent supergiants, the triangles represent giants, and the squares represent dwarfs. Green 
shows the data for the SMC, and red for the LMC. We have distinguished the hottest stars 
(Teff >45,000) by hlled symbols; all of those stars show a signihcant mass discrepancy. 




























